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On the Origin of Cosmic Magnetic Fields 
Russell M. Kulsrud and Ellen G. Zweibel 
Princeton University and the University of 
Wisconsin-Madison 

Abstract 

We review the extensive and controversial literature concerning 
how the cosmic magnetic fields pervading nearly all galaxies and clus- 
ters of galaxies actually got started. Some observational evidence 
supports a hypothesis that the field is already moderately strong at 
the beginning of the life of a galaxy and its disc. One argument in- 
volves the chemical abundance of the light elements Be and B, while a 
second one is based on the detection of strong magnetic fields in very 
young high-red-shift galaxies. 

Since this problem of initial amplification of cosmic magnetic fields 
involves important plasma problems, it is obvious that one must know 
the plasma in which the amplification occurs. Most of this review is 
devoted to this basic problem and, for this, it is necessary to devote 
ourselves to reviewing studies that take place in environments in which 
the plasma properties are most clearly understood. For this reason 
the authors have chosen to restrict themselves almost completely to 
studies of dynamo action in our Galaxy. It is true that one can get a 
much better idea of the grand scope of galactic fields in cxtragalactic 
systems. However, most mature galaxies share the same dilemma as 
ours of overcoming important plasma problems. Since the authors 
are both trained in plasma physics, they may be biased in pursuing 
this approach, but they feel this restriction it is justified by the above 
argument. In addition, they feel they can produce a better review by 
staying close to that which they know best. 

In addition they have chosen not to consider the saturation prob- 
lem of the galactic magnetic field since, if the original dynamo ampli- 
fication fails, the saturation question does not enter. 

It is generally accepted that seed fields, whose strength is of order 
10^^*^ gauss, easily spring up in the era preceding galaxy formation. 
Several mechanisms have been proposed to amplify these seed mag- 
netic fields to a coherent structure with the microgauss strengths of 
the currently observed galactic magnetic fields. 

The standard and most popular mechanism is the alpha-Omega 
mean-field dynamo theory developed by a number of people in the 
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late sixties. This theory and its appHcation to galactic magnetic fields 
is discussed in considerable detail in this review. We point out certain 
difficulties with this theory that make it seem unlikely that this is the 
whole story. The main difficulty with this as the only such amplifica- 
tion mechanism, is rooted in the fact that, on galactic scales, flux is 
constant and is frozen in the interstellar medium. This implies that 
flux must be removed from the galactic discs, as is well recognized by 
the standard theory. 

For our Galaxy this turns out to be a major problem, since unless 
the flux and the interstellar mass are somehow separated, some inter- 
stellar mass must also be removed from the deep galactic gravitational 
well. This is very difficult. It is pointed out that unless the field has 
a substantial field strength, much larger than that of the seed fields, 
this separation can hardly happen. And of course, the alpha-Omega 
dynamo must start from the ultra weak seed field. (It is our philos- 
ophy, expressed in this review, that if an origin theory is unable to 
create the magnetic field in our Galaxy it is essentially incomplete.) 

Thus, it is more reasonable for the first and largest amplification to 
occur before the Galaxy forms and the matter embedded in the field is 
gravitationally trapped. Two such mechanisms are discussed for such 
a pregalactic origin; 1) the fields are generated in the turbulence of 
the protogalaxy and 2) the fields come from giant radio jets. Several 
arguments against a primordial origin are also discussed, as are ways 
around them. 

Our conclusion as to the most likely origin of cosmic magnetic 
fields is that they are first produced at moderate field strengths by 
primordial mechanisms, and then changed and their strength increased 
to their present value and structure by a galactic disc dynamo. The 
primordial mechanisms have not yet been seriously developed as of 
yet, and this preliminary amplification of the magnetic fields is still 
very open. If a convincing case can be made that these primordial 
mechanisms are necessary, more effort will of course be devoted to 
their study. 
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1 Introduction 



It is well established that the universe is filled with magnetic fields of 
very large scale and substantial strength. The fields exist on all scales, 
in planets, stars, galaxies, and clusters of galaxies (Parker 1979). But 
with respect to its origin, the magnetic field in stars and planets is 
secondary, while the field in galaxies is primary. The situation for clus- 
ters of galaxies is not very clear (Carilli 2002); their magnitude and 
structure being rather uncertain. Therefore, the best route to under- 
standing cosmic fields is through discovering their origin in galaxies, 
and in particular in our Galaxy. (Parker 1979, Ruzmaikin et al, 1988, 
Beck et al 1996, Zweibel & Heiles 1997, Kulsrud 1999, Carilli & Taylor 
2002, Widrow 2002, Vallee 2004). 

The idea embraced in this review is: that one has the clearest idea 
of what happens in our Galaxy. If one can not understand the origin 
problem here, then the cosmic origin theory of magnetic fields has to 
be considered incomplete. 

It must be remarked that this choice of reviewing only the work 
on dynamos specifically in our Galaxy, is the choice of the authors 
and represents a certain bias on their part. Generally, reviews of 
galactic magnetic fields discuss the magnetic fields in a great variety of 
extragalactic systems. This in general is justified since, by examining 
the global shapes and properties of fields in external galaxies one can 
form a much better picture of these fields, than by restricting ourselves 
to the field in our Galaxy, in which we see only its more local parts. 
Moreover, the display of these magnetic fields has aesthetic beauty 
which alone should justify this approach. 

However, the authors feel that every one of these extragalactic 
fields represent a very difficult problem from a plasma physics point of 
view. If one wants to understand how all these field in ours and other 
galaxies got started from an extremely weak seed field, one has to first 
deal with fields much weaker than those that can be observed. The 
problem that needs to be overcome is the problem of flux conservation, 
basically a plasma problem. Since the authors are trained plasma 
physicists, they need to know the most basic properties of the plasma 
in which this happened, so these is no better situation to examine 
than our own interstellar medium. Therefore, their choice makes it 
is possible to critically examine the basic plasma physics of galactic 
dynamos in this weak phase and here at home. 

In addition, they do not seriously consider the problem of the sat- 
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uration of the the interstellar magnetic field. In the opinion of the 
authors this problem is really secondary to the origin of the field since 
if the field cannot be amplified by the large amount required to reach 
its present value, the saturation problem does not enter. 

Although it is widely accepted that magnetic fields were not pro- 
duced in the Big Bang, there seems little difficulty with the creation 
of seed fields in the universe, during the period subsequent to recom- 
bination, that is the creation of fields with strengths of order 10"^*^ 
gauss. There are a number of mechanisms that can operate during 
the period of galaxy formation and can generate such fields. The 
Biermann battery (Biermann 1950) is a simple such mechanisrr0. On 
the other hand, the currently observed field strengths are of the order 
of 10^^ gauss. Thus, there is a long way to go between fields with 
these two strengths. Hence, the main problem with the origin of cos- 
mic magnetic fields centers on how the strengths of cosmic magnetic 
fields were raised from weak values of lO"'^'' gauss to the currently 
observed microgauss strengths. 

In discussing this subject of magnetic field origins, we distinguish 
between (1) amplification of fields that are already somewhat strong 
so that the amplification mechanisms can actually by aided by the 
magnetic fields themselves, and (2) amplification of extremely weak 
fields, those whose strength is so weak that they can play no role in 
the amplification mechanisms. The latter are passively amplified by 
mechanisms that are totally unaffected by their presence. 

There is a second division of the problem of amplification that 
concerns the nature of the magnetic field itself. As we will see, it is 
relatively easy to increase the energy of magnetic fields if one allows 
them to be very tangled, changing their direction on very small scales. 
It is much more difficult to produce coherent magnetic fields that 
change their direction only on very large scales, as is the case for the 
magnetic field in our Galaxy. 

Since amplification of the magnetic field energy is relatively easy 
to understand, whether the field is very weak or whether it is strong 
(Batchelor 1950, Kazantsev 1968, Kraichnan & Nagarajan 1967, Kul- 
srud & Anderson 1992, Boldyrev & Cattaneo 2004), the real problem 
of concern for a theory of magnetic field origin is this coherency prob- 

^It is worth noting that magnetogenesis by exotic processes in the early universe has 
also been proposed. Because the nature of these fields by the epoch of galaxy formation 
is highly uncertain, we have ignored them in this review; for good discussions see Grasso 
& Rubenstein (1995) and Widrow (2002) 
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lem, especially for very weak fields. 

Why are we interested in tlic generation and amplification of cos- 
mic fields? There are several reasons. First, until we can be sure 
we understand this problem, we cannot be sure that we really under- 
stand the cosmological evolution of the universe. Second, the actual 
structure of the observed magnetic fields is not very well determined 
because most of the measurements of the magnetic field use tech- 
niques, such as Faraday rotation, that average magnetic fields over 
large distances, (Zweibel & Heiles 1997, Heiles, 1998). If we knew 
theoretically how the magnetic fields were generated, this would give 
extra leverage to determining their local structure. Finally, many of 
the very mechanisms that produce magnetic fields are astronomically 
interesting, and important in themselves. 

Why are galactic magnetic fields astrophysically important? With- 
out their universal presence in the interstellar medium its astrophysi- 
cal properties would be very different. At the present time, magnetic 
fields play a crucial role in the way stars form (Spitzer 1978). They 
also control the origin and confinement of cosmic rays, which in turn 
play important astrophysical roles. Further, magnetic fields are an 
important ingredient in the equilibrium balance of the galactic disc. 

Why is the origin of magnetic fields so difficult to understand? 
First, they are difficult to observe because they are so weak and far 
away. Second, to understand the physics of their origin one needs 
to understand astrophysical plasma physics (Kulsrud 2005, Cowling 
1953), fiuid dynamics, and many other fields of astrophysics. Plasma 
physics on galactic scales is still not a well developed subject and 
the details of how it works are hard to observe. More importantly, 
since the origin of cosmic fields occurs either over the entire life of 
the Galaxy, (Ruzmaikin et al 1988) or possibly in a pregalactic age 
before galactic discs formed, it is very difficult to gain observational 
knowledge of the early generation mechanisms. The early stages in 
other galaxies are observed at large red shift where large distances 
make these observations difficult to make, (Kronberg, 1994, Welter, 
Perry & Kronberg 1984, Watson & Perry 1991, Wolfe, Lanzetta, & 
Oren 1992, Oren & Wolfe 1995.) For these reasons definitive progress 
in uncovering magnetic field origins in the universe has been slow. 

A main goal of this review is to arrive at some sort of conclusion 
as to whether at the stage when the galactic disc formed the magnetic 
field was still extremely weak and the amplification occurred during 
the entire age of the disc, or whether the the fields already were sub- 
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stantial before the galactic disc started to form. If the former is the 
case, we will call the origin galactic and the dynamo generating it 
the galactic dynamo. In the latter case we call the origin pregalactic. 
(We avoid using primordial which suggest a much earlier origin then 
occurs, say, during the protogalactic era.) 

In the next two sections we briefly review the salient observations 
and present a historical introduction to galactic dynamo theory. The 
remainder of the review discusses current theories of magnetic field 
origin and concludes with a critical summary. 

2 Observed properties of galactic mag- 
netic fields 

Our knowledge of galactic magnetic fields rests on four observational 
pillars. Measurements of Faraday rotation and Zeeman splitting give 
the magnetic field perpendicular to the plane of the sky, integrated 
along the line of sight. These effects are direction sensitive, and con- 
tributions from oppositely directed fields tend to cancel each other. 
Observations of the polarized synchrotron continuum, and polarized 
emission and absorption by magnetically aligned dust grains, give the 
line-of-sight integrated magnetic field components in the plane of the 
sky. These diagnostics are sensitive to orientation, not direction, but 
90° swings in orientation also cancel. In addition to line-of-sight av- 
eraging, finite angular resolution of the telescope causes plane-of-sky 
averaging. 

According to these observations, the mean orientation of the mag- 
netic field is parallel to the Galactic plane and nearly azimuthal; the 
deviation is consistent with assignment along the spiral arms (Heiles 
1996). This orientation is consistent with the effects of induction in 
a system with strong differential rotation and a spiral density wave 
pattern (Roberts & Yuan 1970). Assuming that the Galactic halo ro- 
tates more slowly than the disk, induction would act on a vertical field 
so as to produce a reversal in the azimuthal field across the Galactic 
plane (so-called dipole symmetry). The traditional view has been that 
the field does not reverse across the plane (Beck et al 1996), although 
some authors favor asymmetry (Han 2003, Han et al 1997, Han 2001). 
This important question is still open. 

The Galactic field within several kiloparsecs (kpc) of the Sun has 
both mean and random components (Rand & Kulkarni 1989, Han et 
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al. 2006, Beck 2007), with the mean component being of order 1.4-2/y,G 
and the rms field about 5-6/xG. The rms vahie for the random field 
is derived from the assumption that the random field is isotropic and 
from a measurement of its line-of-sight component. There is some 
evidence that the fluctuations are anisotropic, with more power paral- 
lel to the mean field than perpendicular to it (Zweibel 1996, Brown & 
Taylor 2001). There is also evidence that the mean field reverses with 
Galactocentric radius (so-called bisymmetric spiral structure), but the 
locations and frequency of reversal are quite uncertain (see Han & 
Wielebinski 2002 for a review, and Weisberg et al. 2004, Han et al. 
2006, Brown et al 2007 for more recent studies). The discrepant re- 
sults of this important and difficult measurement reflect the high level 
of noise (fluctuations greater than the mean), uncertain distances to 
the background pulsars against which Faraday rotation is measured, 
uncertainty in the location of Galactic spiral arms (Vallee 2005), and 
a possible systematic spatial variation in the field structure which is 
unaccounted for in the models. 

There are some complications in the interpretation introduced by 
the spiral arms perturbing the direction of the magnetic field. Even if 
the unperturbed field is toroidal, these spiral arm perturbations give 
the impression that the global field is aligned along the spiral arms and 
its lines of force are a spiral [Lin, Yuan and Shu, 1969, Manchester, 
1974 (section iv page 642)]. 

How far back in time are galactic magnetic fields detectable? Young 
galaxies and their environment have been probed through the absorp- 
tion lines found in quasar spectra whose redshift is different than that 
of the quasar. These absorption lines are impressed on the quasar light 
as it passes through clouds of gas, and these clouds are interpreted as 
young galaxies. Most of these systems are rather thin and are referred 
to as part of a Lyman alpha forest. However, some of these systems 
are much thicker and display both metallic lines (particularly lines 
of MgH) and very broad damped Lyman alpha lines. The latter are 
referred to as damped Lyman alpha systems. More important for our 
purpose, if the quasar emits polarized radio waves, the plane of po- 
larization of these waves would be Faraday rotated by a significant 
amount if these systems had coherent magnetic fields (Perry 1994). 

To determine this possibility Welter, Kronberg and Perry (1984) 
searched for Faraday rotation in a number of radio emitting quasars, 
and found a definite correlation between those which had a rotation 
measure and those which have metallic absorption lines. The ma- 
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jor difficulty with these observations is the correct subtraction of the 
Faraday rotation produced by the magnetic field of our Galaxy, which 
does not vary smoothly with angular position. Welter et. al. found 
five unambiguous cases. These results were reaffirmed by Watson Sz 
Perry (1991). Since the metallic absorption line systems had red shifts 
that were fairly large (of order two) , these systems probably represent 
young galaxies in an early state of formation. This data was reana- 
lyzed by Wolfe and his collaborators (Oren and Wolfe, 1995), but in 
a different manner. They restricted their analysis to 61 quasars with 
Mgll absorption lines and separated out, as a subclass, 11 of these that 
also had damped Lyman alpha lines. They found that the incidence 
of Faraday rotation in the 11 damped systems was higher than that in 
the remaining 50 undamped systems with a 99.8 per cent confidence 
level. 

In this analysis they concluded that the errors introduced by the 
Faraday rotation in our Galaxy were larger than those assumed by 
Welter et al by a factor of three. Thus, they only found two cases in 
which they were certain that there was Faraday rotation. 

The detection of only two cases with definite intrinsic rotation 
measures seems to make a weak case for extragalactic fields in these 
damped systems. But these cases were those systems with the lowest 
red shift. For the other systems of larger redshifts z, any intrinsic 
Faraday rotation produced by them is diluted by a factor of (1 + z)^. 
This is because the frequency of the radio waves passing through them 
is higher by the factor (l + z) and the amount of Faraday rotation de- 
creases with the frequency squared. Thus, the other members of the 
damped class could very well have magnetic fields of the same strength 
as the lower red shift members without producing a detectable signal. 
This bolsters their correct subtraction of the galactic component of 
the Faraday rotation. Taking this into account, Oren and Wolfe con- 
clude with 98 per cent confidence that all such systems have Faraday 
rotation and substantial magnetic fields. 

Another important window on the history of the Galactic magnetic 
fields over cosmic time is provided by analyzing the chemical composi- 
tion of the atmospheres of the oldest stars in the Galactic halo (Zweibcl 
2003). As a result of observations from the Hubble space telescope the 
chemical abundances of these very early stars have now been analyzed. 
The light elements lithium, beryllium, and boron have been found in 
even the oldest stars, e.g. those with 10"^ times solar abundances. In 
addition, the amount of beryllium and boron in them increases with 
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the iron abundance and in fact is directly proportional to it. Since 
the early stars arc produced from the interstellar medium their com- 
position should reflect that of the interstellar medium. (Primas et al. 
1999, Duncan et al. 1998, Garcia-Lopez et al. 1998, 1999) 

Now, it is known that no beryllium was produced in Big Bang 
nucleosynthesis, and that it is very difficult to make it in stars, since it 
quickly burns up. The leading theory is that it was made by cosmic ray 
nucleosynthesis, (Meneguzzi, Audouze &; Reeves, 1971, Reeves, 1994, 
2007). The situation for boron is ambiguous, because this element can 
also be produced by neutrinos during Type II supernova explosions; 
see Ramaty et al. 1997. We will still include boron in the arguments 
for cosmic rays and magnetic fields early in the history of the Galaxy, 
but this caveat should be kept in mind. 

The linearity of the Be and B abundances with iron is explained 
by assuming their creation is by spallation of the low energy (tens of 
MeV) carbon and oxygen cosmic rays. If it were due to the comple- 
mentary process, spallation breakup of interstellar carbon and oxygen 
nuclei, one must take into account that the latter themselves are pro- 
duced by stellar nucleosynthesis and supernovae. and their abundance 
should increase proportionally to the amount of iron produced in su- 
pernovae. Thus, their abundance should increase with that of iron in 
the interstellar medium, and the abundance of these light elements 
should be quadratic with iron abundance or in time. 

However, to get linearity in time with spallation of carbon and oxy- 
gen cosmic rays, one needs to assume that the composition of cosmic 
rays must not change with time. This would seem to be a stumbling 
block since cosmic rays are assumed to be accelerated by shocks from 
interstellar-medium nuclei. Therefore, they would be expected to also 
reflect the changing abundances in the interstellar medium and would 
be expected to increase their abundance in carbon and oxygen with 
time. 

Ramaty and others (Ramaty, Lingenfelter, Kozlovsky (2000), Ra- 
maty, Scully, Lingenfelter, 2000a, Lingenfelter, Higdon, Ramaty, 2000b, 

Ramaty, Tatischcff et al. 2000) argue that the acceleration of cosmic 
rays occurs mainly inside superbubbles, and that the material inside 
these superbubbles is the material that has just emerged from the 
supernova generating the superbubble. This freshly produced matter 
has not been diluted with the rest of interstellar medium in the super- 
bubble region where cosmic ray acceleration occurs. Thus, the relative 
abundance of different cosmic ray nuclei should be constant in time 
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and determined by the undiluted material emerging from supernovae. 
On the other hand, iron in the interstellar material from which stars 
emerge has been diluted, since it was produced in supernovae. There- 
fore, its abundance relative its solar abundance increases with time at 
a constant rate determined by the rate of supernova explosions. 

This constancy of the chemical abundance of cosmic rays has been 
supported by detailed numerical simulations of Ramaty and others, 
and does lead to an explanation of the observations that the abundance 
of beryllium and boron is directly proportional to the abundance of 
iron in old stars (Garcia-Lopez et al. 1998, 1999) 

Now to explain the numerical value of this ratio, the cosmic ray 
intensity at tens to hundreds of MeV per nucleon in the early Galaxy 
had to be the same then as now. Zweibel (2003) has shown that 
magnetic fields several orders of magnitude weaker than now, suffice 
for cosmic ray acceleration and diffusive propagation at these energies. 
On the other hand, if the field is very weak such a cosmic ray intensity 
might not be confined by a very weak field. This is because, if most 
of the mass in the interstellar medium was in the form of discrete, 
cold clouds, as appears to be the case today, then a high cosmic ray 
pressure between the clouds, which must be confined by magnetic 
tension, would inflate the field lines to infinity (Parker 1979). 

This can be quantified by a simple two dimensional model. Let 
the cloud distribution be two dimensional and have a scale height H 
and the clouds be on planes separated by a mean distance £. Let 
the cosmic ray pressure be proportional to the magnetic field strength 
squared by a factor f3/a. Then the model shows that the lines bow 
out above their average height in the clouds by a factor 

1 

cos ^y{P/a + l)£/H 

Thus, if the field is very weak, P/a is very large and there is no 
solution, implying the field lines would bow out to infinity. 

Although the amount of cold interstellar material in the early 
Galaxy was doubtless lower than today because of the lower metal- 
licity, thereby somewhat easing the cosmic ray confinement problem, 
we are faced with a situation in which primeval galaxies already must 
have had substantial magnetic fields, at a stage too early to have been 
produced by a conventional dynamo. 

Finally, the problem of the abundance of the Li^ isotope which is 
not at all linear with that of iron is entirely up in the air. (Reeves 
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2007, Ramaty, Tatischeff et al. 2000). Li^ also is not produced by Big 
Bang nucleosynthesis and, because it depletes so easily, can only be 
produced with great difficulty in stars. Thus, as long as we do not 
understand the process that produces lithium 6, we cannot be sure 
that this process does not in some way produce beryllium and boron 
as well, without the aid of cosmic rays. Until the Lithium problem 
is resolved, we cannot be certain that the beryllium boron argument 
proves that the origin of Galactic field is pregalactic. 

It is interesting that boron has been observed in other galaxies at 
large red shifts (Prochaska et al 2003). Thus, if the ideas concerning 
the origin of boron hold up, this provides even stronger evidence that 
magnetic fields are already present at the formation time of galactic 
discs and their origin is pregalactic. 

3 Summary of our present understand- 
ing of cosmic magnetic field origins 

Origin theories divide into two further parts: (1) the origin in small 
bodies such as planets and (2) the origin in larger bodies such as 
galaxies. This division occurs because, for a large body such as a 
galaxy, the resistive decay lifetime of its magnetic fields is much longer 
than the age of the universe. For these bodies a so-called fast dynamo 
is required. On the other hand in small bodies the decay lifetime is 
much shorter and the required dynamo is called a slow dynamo. For 
the Earth the problem of where the field came from, and how it is 
sustained against resistive decay is easier than the same problem for 
fields in large bodies (Parker 1979, Spitzer 1978). If a magnetic field 
starts to decay, the inductance of the body in which it resides produces 
a backvoltage (or E.M.F.) that keeps its currents flowing against its 
resistance. Thus, one can roughly say that the lifetime of the magnetic 
field against decay (if no other mechanisms are present) is its total 
inductance divided by its total resistance. Inductance is proportional 
to the size of a body, L while resistance is inversely proportional to 
L so that the lifetime is proportional to L^. Thus, the Earth has a 
relatively short magnetic decay time of order a few tens of thousands 
of years, while that of the Galactic disc is many orders of magnitude 
longer than the age of the universe. Thus, dynamo mechanisms aside, 
the Earth's field would decay away in a time very short compared to 
its age, and therefore, there must be a mechanism to sustain it, just 
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as to sustain a current in a laboratory circuit, one needs a battery or 
dynamo (Parker 1955, Elsasscr 1946, 1950.) 

On the other hand, while one need not worry about the galactic 
field decaying because of its enormous inductance (Fermi 1949), one 
has to worry about how to get the currents started in the interstellar 
medium to produce the galactic field: as the currents rise the back 
voltage is so large that a very strong generator is needed to balance 
the back voltage. This turns out to be the essence of the problem 
behind the origin of galactic fields, (Hoyle, 1958). 

This discussion of magnetic field generation and decay is not quite 
correct. It treats the bodies as being static and at rest, while in both 
cases the bodies are either fluid or gaseous with motions generated 
by their dynamics. When a fluid moves across a magnetic field, an 
electric field E = — v x B/c exists in the frame in which the veloc- 
ity is measured, and this electric field results in the dynamo action 
that is needed to balance the magnetic field against resistive decay in 
case of the Earth, and balance inductance in the case of the Galaxy. 
The problem is to find a reasonable fluid velocity that would properly 
balance the inductive and resistive effects that must occur during the 
evolution of the field (its decay or growth). 

In 1955 Parker was able to find such velocities and to propose 
a model containing, non axisymmetric motions, that explained how 
the Earth's field could be sustained against decay (Parker 1955). To 
do this, he built on the work of many others (Elsasser, 1946, 1950). 
Parker found that a dynamo should exist in the liquid core of the Earth 
and the fluid motions producing this dynamo action arc a combination 
of differential rotation of the core, and a multitude of rising and falling 
convective cells, that are twisted by the Coriolis force of the Earth's 
rotation. 

His solution was gradually improved in the next decade (Backus, 
1958), until finally, in 1966, Steenbeck, Krause, and Radler (1966) 
developed a refined theory for dynamos, the mean-field theory, which 
consist of such turbulent motions 

Once this theory was accepted as correct, it was applied to the 
problem of the origin of the Galactic magnetic field. Parker in the 
United States (Parker, 1971a) and Vainshtein and Ruzmaikin in Rus- 
sia (Vainshtein &; Ruzmaikin, 1972) in the early seventies showed that 
motions, similar to the terrestrial motions, exist in the Galaxy, and 
that they could overcome the inductance problem and generate the 
Galactic field. 
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Although the dynamo mechanism appears viable, it still docs not 
resolve the primary origin question: was the galactic field actually 
produced, from an initially weak field, during the galactic lifetime, or 
was it generated earlier, before the galaxy formed? That is, for the 
latter case, was a fairly strong field already present at the onset of 
formation of the galactic disc, and then, its final form determined by 
this dynamo action? (See Kulsrud, 1999). 

In the last analysis, this question can only be answered by observa- 
tions. As pointed out above, direct observations of the early magnetic 
fields are very difficult. At the present moment there are two ways to 
make progress on the origin problem: First, a careful theoretical anal- 
ysis of the proposed dynamo model and its assumed velocities; second, 
an analysis of the effects of an early magnetic field on independent as- 
trophysical problems such as star formation and cosmic ray properties, 
discussed in §2. If the field changed drastically from being very weak 
to its present value, these properties, which are subject to direct ob- 
servations, should also change. The stellar luminosity function would 
be expected to change from its early form to its present form. Also, 
the properties of cosmic rays should drastically change. Their energy 
density should either increase from an initial small value to its present 
value, or should remain constant. Its variation can be inferred from 
the measurement of the abundances of the isotopes of elements that 
can only be produced by cosmic ray spallation, as discussed in the 
previous section. 

In the seventies, when the galactic dynamo origin of magnetic fields 
was first proposed, the velocities needed for the dynamo were not well 
known and the estimates made of their magnitudes led to too weak a 
dynamo to be able to amplify the strength of seed fields to the current 
strength magnetic field, over the finite age of the Galaxy. More precise 
estimates were made in a series of papers by Ferriere in the nineties 
(Ferriere, 1992a,b 1993a,b, 1996, 1998, Ferriere and Schmitt, 2000). 
In these papers she pointed out the importance of the expansional 
motions of random super bubbles; [giant bubbles produced by multiple 
supernovae that arise in stellar clusters (McCray and Snow 1979)]. 

These motions last long enough, and are sufficiently fast for the 
necessary amplification of the fields. However, they work only if signif- 
icant magnetic flux can be expelled from the Galaxy. This is possible 
only if the magnetic field lines can separate from the heavy plasma and 
evaporate away leaving the plasma in the superbubble to fall back into 
the Galactic disc. But arguments that this is possible assume that the 
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magnetic field is already fairly strong (Brandenburg et al, 1995, Moss 
et al 1998). Thus, the field strengths they require would be inadequate 
to make the dynamo work in the early stages of amplification if the 
initial field at the beginning of the Galaxy were a weak seed field. 

To be sure the dynamo could work during early galactic times, one 
needs to consider the possibility that a field of considerable strength 
and of pregalactic origin, existed prior to the galactic disc. If the 
pregalactic field is strong enough to solve the problems associated 
with the galactic dynamo theory, then one may have a double origin. 
First, a pregalactic magnetic field of moderate strength must develop, 
and then a subsequent amplification of it to its present strength and 
configuration be produced by the conventional galactic dynamo. 

This pregalactic field could arise as a natural part of the galac- 
tic formation process. For example, it is likely that the formation 
of the first massive stars preceded the formation of galactic disks, 
occurring instead in the ~ IO^Mq halos that appear to be the first 
bound structures formed (Tegmark et al 1997, Gilmore et al. 2006). 
Any magnetic field present in the ambient gas or in the supernova 
ejecta would have been amplified by turbulence in the halo gas. Thus, 
the disks that formed as the halos merged, could have had moder- 
ately strong magnetic fields from the outset. Turbulent amplification 
can also occur during the gravitational formation of larger structures 
where very strong large scale turbulent velocities are present (Ryu et 
al 1993, Kang et al. 1994). 

These velocities are sufficiently large that sufficiently strong mag- 
netic fields can be generated to provide a pregalactic origin (Kulsrud 
& Anderson 1992). If these fields are formed on the scale of the largest 
turbulent eddies which drive them, this would provide the coherence 
needed for a pregalactic origin of the galactic field. However, the the- 
ory of these random turbulent fields is still under development and a 
number of tricky plasma problems have arisen involving them. (See 
§14.) Simulations of turbulent amplification show a pileup of magnetic 
power at the resistive scale in situations where the resistive scale is 
much smaller than the viscous scale (Schekochihin & Cowley, 2004). 
Even smaller scale fluctuations may arise in the coUisionlcss case, be- 
cause of anisotropic plasma pressure (Schekochihin et al 2005b). It 
appears that the fate of a primordial theory based on the protogalac- 
tic physics may hinge on the resolutions of these plasma problems. 
The question of generation of significant fields at earlier epochs in the 
universe than that of galactic disc formation are still being debated. 



15 



It has also been proposed that the pregalactic fields could arise 
from the fields observed to be present in giant radio jets, (Daly &: 
Loeb, 1990, Colgate & Li 2000, Kronberg et al. 2001, Furlanetto & 
Loeb, 2001). So far, none of these proposals has been developed to a 
sufficient extent to judge whether they lead to a viable solution. 



4 Basic equations for magnetic field 
evolution 

Nearly all theories of magnetic field origin rely on the magnetohy- 
drodynamics (MHD) description of plasma physics. It turns out that 
the larger the system the better these equations describe magnetic 
field evolution. These equations essentially parallel the ordinary elec- 
tromagnetic equations. One equation describes how the motion of a 
plasma affects the magnetic field, and the other describes how the 
magnetic field affects the motion of the plasma. Although these equa- 
tions are familiar, we sketch their derivation, since this derivation 
demonstrates the accuracy of these equations for our purpose. 

We start first with the effect the plasma motion has on the mag- 
netic field (Cowling 1953). This effect results solely from the fact that 
in the frame of the moving plasma the electric field E', to all intents 
and purposes, must balance the resistive term, so that we have 

E' = E+^ = r?j, (1) 

c 

where v is the plasma motion, the left hand side is the electric field 
in the frame of the plasma, rj is the resistivity and j is the plasma 
current density. The right hand side is essentially zero since from 
Ampere's law the current density is inversely proportional to 1/L, 
and on galactic scales is usually very small, while the resistivity is the 
same as that in a laboratory plasma. Combining this with Faraday's 
law dB/dt = — cV X E, and Ampere's law V x B = 47rj/c, we get 

f = Vx(vxB) + |^V^B. (2) 

We drop the displacement term from Ampere's law and, for simplicity 
take 77 to be a constant, to get the last term. 

The neglect of the resistivity term in galactic physics becomes ob- 
vious if we approximate V^B as B/L^ and set the last term equal to 
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^/Tdecay For a plasma at a temperature of 10^ degrees Kelvin, the 
resistive term rjc^/An ^ lO^cm^/sec so 

Tdecay ~ lO'^L^, years , (3) 

where Lpg is the scale size in parsecs. Since the galactic scale is many 
parsecs, the incredible smallness of the resistive term in equation (2) 
is evident. Even if the magnetic field is extremely tangled, so that its 
length scale is reduced to a fraction of a parscc, the resistivity term 
is still very small. (For L > 10^^ cm the effective decay time is still 
longer than the current Hubble time.) 

Thus, in almost all cases, it is permissible to drop this term, and 
this is valid to an extraordinary degree. There are exceptions to this 
in cases where extremely thin current layers form. Also there are 
additional, very weak terms in Ohm's law that we come to later, which 
can play a role in producing seed fields. An example of these terms is 
the so called Biermann battery. 

For the bulk of the theories of magnetic field origin the ideal mag- 
netic differential equation 

^ = V X (V X B) (4) 

is sufficiently correct. 

The ideal equation has an important implication for cosmic mag- 
netic field origin theories. It implies that magnetic field lines have a 
physical reality beyond their mathematical description for a magnetic 
field. (Their mathematical definition is: that their signed direction 
gives the signed direction of the magnetic field and their density can 
be assigned so as to give the magnetic field strength.) 

Their reality consists of the fact that the magnetic field lines cannot 
be created or destroyed once they are embedded in the large scale 
plasma. They are bodily carried by the motion of the plasma in which 
they are embedded, and continue to properly represent the magnetic 
field. The argument that flux lines have a physical reality is standard, 
and is called flux freezing. It is discussed in many textbook on plasma 
physics (see, for example, Newcomb 1958, Kulsrud 1964, Alfven Sz 
Falthammer, 1963, Moffatt, 1978, Parker 1979, Kulsrud 2005). 

The mechanism for the creation and destruction of fleld lines is the 
resistivity, which is negligibly small. But, by its very definition, an 
origin theory starts off with a very weak field and produces a strong 
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field. Because of flux freezing, this has to be done without changing 
the number of lines of force. 

It is a fact that stretching a tube of force by motions that preserve 
its volume, does increase the local field strength. This stretching oc- 
curs in a finite volume by doubling the line back and forth as in Figure 
3 of §14, so that the mean value of the field strength, taken over the 
volume, does not change and, thus, does not lead to an increase in the 
coherent field. 

It turns out that standard dynamo theories actually work this way. 

They fold lines of force back and forth increasing the pointwise field 
strength. The dynamo works in such a way that the field lines pointing 
in one direction are near the center of the plasma volume and the 
field pointing in the reverse direction are concentrated on the outside. 
Thus, the magnetic field becomes coherent over roughly the inner 
half of the volume with flux of one sign, and coherent over the outer 
half with flux of the other sign. In §6, we show how this remarkable 
process happens when we discuss the conventional a—U theory for the 
exponential growth of a magnetic field in the galactic disc. However, 
to actually produce a strong magnetic field of constant sign over an 
entire volume, the magnetic flux in the outer region must be expelled. 

A leading question then becomes whether the flux can actually be 
expelled to complete the galactic dynamo action. If one were content 
with a field of reversing sign, constant over one half the volume, but 
with zero total net flux, one could avoid this expulsion problem. How- 
ever, the observed magnetic field in our Galaxy appears to have net 
flux, so that such a solution to the origin problem for our Galaxy is 
not available. The same difficulty occurs in other galaxies as well. 

Moreover, there is an issue of scale. The typical strength of seed 
fields at the onset of the disc is 12 - 14 orders of magnitude below 
present day galactic fields. In order to bring them to their present 
strength, the field lines have to be stretched by 12 - 14 orders of 
magnitude as well. This would seem to imply that the field lines 
fold down to length scales that are also smaller by 12 - 14 orders 
of magnitude or, without any flux expulsion, down to the resistive 
scale. But by invoking flux expulsion at every doubling of the flux 
one can avoid this problem. In fact, only tiny pieces, of held lines 
~ 10~^^ — 10~^^ of the total initial line length are stretched by this 
large order of magnitude while the rest of the pieces of the line are 
removed from the disc. 

The other half of the MHD picture of how magnetic fields affect 
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plasma motions is of less interest in the origin problem since the mag- 
netic fields we consider are too weak to affect the plasma motions. 
Their effect is the introduction of a magnetic force term 



into the hydrodynamic equation. As remarked in §1, we do not con- 
sider the saturation problem in this review since it clearly does not 
arise unless a successful build up of the field to large strength is not 
achieved. 



5 Cowling's theorem and Parker's dy- 
namo 

We have remarked that the motion of a plasma acts like a dynamo 
in that, as a plasma moves perpendicular to a magnetic field, a sub- 
stantial electric field, — v x B/c, arises in the laboratory frame. This 
was first pointed out by Larmor (1919). In fact, this is just how a 
dynamo in a power station generates electricity. Larmor expressed 
the hope that velocities could be found which could generate astro- 
physical magnetic fields. The first attempts to create magnetic fields 
in this way in spheroidal bodies naturally assumed that the velocities 
and fields were axisymmetric. 

Axisymmetric fields are generally broken into poloidal components, 
whose lines of force lie in meridional planes, and toroidal components, 
whose lines of force form circles about the axis. Imagine starting with 
a pure poloidal field with a zero toroidal component, and imagine that 
the plasma is differentially rotating with a rotation rate that varies 
along a poloidal field line. It is clear that this rotation will drag 
different parts of the line in the toroidal direction at different rates, so 
the line will not stay poloidal but will develop a toroidal component. If 
this continues, the toroidal component will grow stronger and stronger 
at a linear rate in time. 

Prom this it is easy to see that a toroidal magnetic field is easy 
to generate from a poloidal field by differential rotation, (see Kulsrud 
2005). On the other hand, if one starts with a pure toroidal field, it 
is impossible to produce a poloidal component. Rotation will have no 
effect on the field and any symmetric poloidal motion will just move 
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the toroidal field lines around leaving them as circles, although the 
circles will change their position. 

Cowling (1934) showed that this is the general situation: if one 
starts with a magnetic field that has both axisymmetric toroidal and 
poloidal components, then axisymmetric motions can increase the 
toroidal flux, but they will just leave the total poloidal flux unchanged. 
This result is known as Cowling's theorem and basically defeats any 
attempt to exponentially amplify a very weak field to a strong one 
by axisjnumetric motions. (The original purpose of the theorem was 
show it was impossible sustain the Earth's field against its relatively 
rapid resistive decay. However, it applies equally to any attempt to 
produce a galactic dynamo that starts with a weak axisymmetric co- 
herent field, and tries to increase it to a coherent field of order 10~^ 
gauss by axisymmetric motions.) To increase the toroidal magnetic 
field linearly by purely differential rotation would take 10^^ differen- 
tial turns, while the Galaxy has only rotated about fifty times. (It 
should be remarked that it is possible to amplify a non axisymmetric 
magnetic field by symmetric velocities, but this does not seem to have 
been attempted for a galactic dynamo.) 

Thus, to exponentially increase the flux of an initial small seed field 
to the present Galactic value, or to sustain the Earth's field against 
decay, it is necessary to deal with nonaxisymmctric motions. This 
would seem to require complicated numerical simulations, which were 
beyond the power of computers of the 1950's. 

However, as previously mentioned, Parker (1955) did come up with 
a physically correct model containing non axisymmetric motions that 
could successfully, either amplify the Earth's field, or sustain it against 
resistive decay. He considered random small convection cells of rising 
and falling fluid in the Earth's liquid core. Imagine, as above, that 
there is a purely toroidal magnetic field inside the Earth's core. Then 
it is true that by itself a rising convection cell would distort this field, 
but only in a plane not producing any net poloidal component. How- 
ever, due to the Coriolis force of the Earth's rotation the distortion 
of the toroidal field would be twisted into the poloidal plane and a 
poloidal component would emerge. If one only had this convection, 
the poloidal field would not grow significantly. However, when com- 
bined with the differential rotation of the Earth's core, the poloidal 
distortion would produce a change in the initial toroidal component 
of the field. (This differential rotation in the Earth's core also results 
from Coriolis forces.) 
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It turns out that the change in the toroidal field is actually an 
increase. Thus, one has the possibility of exponential growth. In the 
lifetime of a convection cell the poloidal component increases by an 
amount proportional to the toroidal component, and the toroidal com- 
ponent increases by an amount proportional to the poloidal compo- 
nent. Since these two increases are coherent, both components grow 
exponentially. Parker showed that all the signs are consistent with 
growth. The downward motions produce changes with the same sign 
as the upward motions, and the same changes apply to both the north- 
ern and southern hemispheres. Thus, all of the random convection 
cells act to increase the mean magnetic field of the earth in the same 
direction. 

In detail, a rising convective cell rotates in the opposite direction to 
the Earth due to the Coriolis force. This is because it laterally expands 
increasing its moment of inertia. This leads to a twist in a distorted 
toroidal field line with the upper part of the line gaining a northward 
component and the lower part a southward component. The Earth 
rotates slower at larger distance from its axis, so the differential ro- 
tation acting on this new poloidal field loop produces an additional 
toroidal field in the same direction as the original field reinforcing it. 
A sinking convection cell contracts and rotates oppositely, yielding a 
loop with the same sense as the rising cell, and also reinforcing the 
toroidal component. 

The cells in the southern hemisphere act the same, reinforcing 
the toroidal field in the southern hemisphere which is in the opposite 
direction to that in the northern hemisphere. Thus, because the sign of 
the Coriolis force is properly correlated with the differential rotation 
of the Earth's core, the convective cells all lead to a growth in the 
toroidal field, or compensate for the resistive decay and lead to a non 
decaying terrestrial magnetic field. 

This theory was intuitively satisfactory, and was also subsequently 
supported by the more rigorous analysis of Backus (1958), and by later 
numerical simulations (Glatzmaier k, Roberts 1995, Ogden et al 2006). 
It can most easily be illustrated, in the galactic case, by the example 
of a supernova exploding in a rotating galaxy. See figure 1, (Kulsrud 
2000) . The convection cell is illustrated by the exploding remnant that 
is forced to rotate backward by a Coriolis force induced by Galactic 
rotation. The role of differential rotation is also supplied by Galactic 
rotation. In this Galactic case, the lines need to be expelled to infinity, 
in contrast to the terrestrial case where they need merely be expelled 
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into the nonconducting mantle and effectively disappear. 

However, to successfully apply Parker's theory to the Galactic mag- 
netic field problem a more systematic treatment of turbulent amplifi- 
cation of magnetic fields is required. This was provided by the mean 
field dynamo theory of Steenbeck, Krause and Radler (1966). They 
treat the turbulent motions in a more general way than Parker treated 
the Earth, replacing the convective cells by a quantity they called ki- 
netic helicity, v • V x v. They denoted its effect on magnetic fields by 
a where 

a = -- < V • V X V >, (6) 

where the angular brackets denote turbulent ensemble averages and r 
is the decorrclation time of the turbulent motions. (To obtain a scalar 
a, isotropic turbulence must be assumed.) Note that helicity has the 
same sign of twist relative to the motion v as do Parker's convective 
cells, and serves the same function in amplifying the field on which 
they act. (rV x v) corresponds to the angular twist of the convective 
cell during a single decorrclation time. 

Because the normal resistivity is usually very small, they also in- 
troduce a turbulent resistivity which they called /3 where 

/? = ^<v2>. (7) 

This turbulent resistivity does not act the same as normal resistivity in 
that it produces no dissipation or change in field line connectivity. It is 
an effective mixing term produced by random turbTilent motions, and 
smooths out the fields on a larger scale than the turbulent motions. 
It can be justified as follows: r v is the displacement Ar, so that 
/? K, {i^vf/r is the expression for a random walk. 

6 The alpha-Omega disc dynamo 

The mean field dynamo theory introduced by Steenbeck et al (1966) 
was a very important step towards applying dynamo action in many 
contexts: the Earth, the Sun, stars, galaxies, and clusters of galaxies. 
We present it in the context of our Galactic disc, which is very thin. 
This allows us to treat the dynamo locally in radius and angle and 
makes it one dimensional with the significant space coordinate perpen- 
dicular to the disc. This is consistent with the general Galactic field, 
and is the case most relevant to the Galactic dynamo (Ruzmaikin, 
Shukurov, & Sokoloff 1988). 
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We start with a very weak field such that the motions can be con- 
sidered to be independent of the magnetic field, a limit called the 
kinematic limit. Thus, one is only concerned with the magnetic differ- 
ential equation (2). Because resistivity is so small, we can start with 
the ideal equation (4). In this equation we split the velocity v into 
two parts, a part associated with the turbulent random motions Sv, 
and a part U which is smooth and coherent. In the Galactic case U 
is the differential rotation of the Galaxy. Thus, 

V = U + 5v. (8) 

The magnetic field is also broken up into its mean part B and its 
random part 5B, so 

B = B + (5B. (9) 

Substituting these into equation (4) and ensemble averaging over the 
turbulence we get, 

^ = Vx(UxB)+Vx(((5vX(5B)). (10) 

Thus, an extra term associated with the random velocities has been 
added to the magnetic differential equation for the smooth field. This 
term can break the ideal flux constraint for the mean field B, but not 
for the true field B. 

To complete the process of obtaining a dynamo equation, one has 
to solve a magnetic differential equation for (5B in terms of S\. This 
is carried out by what is known as a quasilinear expansion (Sagdeev 
and Galeev 1969), which is presented in many places and will not be 
given here [see Ruzmaikin, Shukurov, & Sokoloff, 1988; Parker, 1979; 
Krause & Radler, 1980; Moffatt 1978; and Kulsrud 2005). 

The result for the mean induction generated by the fluctuating 
fields can be written 

{SvxSB) =aB-pV xB. (11) 

[In general, a and /? are tensors, but here we understand them to be 
scalars and defined through equations (6) and (7)]. Substituting this 
result into the above equation gives the famous mean field dynamo 
equation 

BB — — — 

— = Vx(UxB)+Vx (aB) + pV^B. (12) 
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Applying this equation to a local part of the Galactic disc, in 
which we introduce cylindrical coordinates r,9, z, and keeping only 
derivatives in the thin vertical direction, z, we get 



dBr 




d'^Br 



dt 

dBe 



dt 



= -nBr + p^-^. 



(13) 



(We have dropped the bars; the magnetic field B is understood to be 
only the mean field.) 

Wc have substituted the Galactic rotation ri}9 for U, and used the 
fact that rO, is a constant to find its derivative dQ,/dr = — O/r. We 
have also dropped the term dz{aBr) in the equation for Bg, which is 
small relative to the differential rotation term. 

The usual procedure is to look for a growing mode, proportional 
to e^* and to solve equations (13) as an eigenvalue problem. The 
boundary conditions that are customarily invoked are that the disc 
is confined to a region —h<z<h. On the presumption that the 
diffusion coefficient, (3, is very large outside of the disc, the magnetic 
field is taken to be zero foi \z\ > h, and therefore the boundary condi- 
tions are: Bj. and Bq are zero at \z\ = h. These conditions are termed 
the vacuum boundary conditions. As we will show later (in §9), the 
assumption of these conditions is the critically important point in the 
disc dynamo origin for the Galactic field. As discussed in §3, in or- 
der to preserve flux conservation, flux of one sign must be removed 
from the disc. This is exactly what these boundary conditions accom- 
plish. However, the weak point in the Galactic dynamo theory is the 
assumption that /3 is large outside the disc. This too is discussed in 



Although the assumption that P is constant is non objectionable, 
a cannot be constant since it has to be of opposite sign above and 
below the disc. This can be seen physically if one interprets a as the 
effect of twisting convection cells. 

In fact, consider a cell rising above the midplane of the disc, that 
is z and Vz > 0. As it rises it expands, and the Coriolis force of galac- 
tic rotation rotates it backward relative to normal Galactic rotation. 
(This follow easily if one considers that the moment of inertia of the 
convective cell is increasing so its angular rotation relative to a flxed 
frame must decrease to preserve its angular momentum.) The same 
backward rotation applies for a descending cell below the galactic mid- 



9. 
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plane. However, for this cell Vz < 0, and a ~ — v - (V x v), the hcUcity, 
reverses sign. Inspection of equation (6) of the last section shows that, 
since Galactic rotation is positive around z, helicity is negative above 
the midplane and a is positive. 

Given the uncertainty in the properties of interstellar turbulence, 
it is difficult to give an explicit form for a, so various forms have 
been tried in order to explicitly solve the eigenvalue equations. The 
simplest form a = aoz/h, is probably as good as any. 

To solve the eigenvalue equations (13), they are cast into dimen- 
sionless form. 

z' = z/h, t' = pt/h^, y = 7/iV/3, 
Be = B'ip/hao), Br = B',iP/hao). (14) 

so that the dimensionless eigenvalue equations become 

I'B' = -Dft + ^ (15) 

where the dimensionless number 

^ = -^. (16) 

determines the eigenvalue 7'. (Since Q, h, ao and /3 are positive, D is 
negative.) 

D is known as the dynamo number. It represents competition 
between the growth terms (shear and the a effect), and diffusion (the 
P effect) and determines whether 7 is real and positive. There is a 
negative critical value Dc such that if D < Dc, then 7 is positive and 
the magnetic field will grow exponentially. The dimensional growth 
time, if 7' = 1, is h'^/p. 

Following Parker or Ruzmaikin, Shukurov and Sokoloff, one ob- 
tains an approximate idea of the size of /? in the interstellar medium, 
by taking the random turbulent velocity Sv « lOkm/sec , and the cor- 
relation length 5vT ~ 100 pc. This gives /3 ~ 1.5 x 10^^ cm^/sec. A 
possible value for h is 300 pc. This gives for the dimensionless growth 
time 1.5 x lO^^sec = 5 x 10^ years. That is to say, if the eigenvalue 
problem leads to 7' = 1, an initial seed field will exponentiate every 
500 million years, and there could be roughly twenty efoldings during 
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the age of the Galactic disc. This would raise an initial field by a 
factor of 10^ from say a field of 10~^^ gauss to 10~^ gauss. 

It should be noted that in the Galaxy today, while 300 pc is proba- 
bly a reasonable scale height for a, the observed magnetic scale height 
appears to be 1500 pc. However, since the field is anchored in the 
thinner interstellar mass, 300 pc is the better value. The field at 1500 
pc. is confined by magnetic tension which has not been taken into 
account in the simplified dynamo theories. Using a larger value of 
h would decrease the rate of diffusion, increase \D\, and increase the 
growth rate of dynamo waves. 

Since the a — dynamo leads to exponential growth of the mean 
field, it must eventually saturate. This is expected to happen by 
nonlinear quenching of a and /3. Many papers have been devoted to 
exploring this quenching. (For example, see the excellent papers of 
Blackman & Field, 1999, 2000 or Kleeorin et al. 200, 2002, 2003). 
The crux of the issue is whether the dynamo is quenched when the 
mean field B comes to equipartition with the turbulence (which is ap- 
proximately the situation in our Galaxy) , or whether small scale fields 
suppress the dynamo long before B reaches equipartition. We return 
to the issue of the small scale fields in §14. For an extensive review 
the saturation problem consult the Physics Report of Brandenburg 
and Subramanian (2005), and also Shukurov (2004). 

Proceeding with the mean field analysis and examining equations 
(13), we see that they are symmetric in z so that there are solutions 
with two different parities. In one both magnetic field components are 
even in z, and in the other they are both odd. The first solution is 
called the quadrupole solution. It has a much smaller critical dynamo 
number than the odd one which is called the dipole solution. The 
quadrupole solution is the one whose parity agrees most closely with 
the present Galactic field, which appears to be symmetric with respect 
to the midplane. However, this is the solution in which the net flux in 
the disc changes and, therefore, requires the flux expulsion provided 
by the vacuum boundary conditions. Its critical dynamo number is 
about —13 for the form wc have chosen for a{z). 

If one had closed boundary conditions with 
dBr/dz = dBg/dz = 0, no such flux problem arises, and one would 
not need flux expulsion. For this case, the critical number for the 
dipole solution is —4, and a growing quadrupole solution does not 
exist. Thus, knowing nothing about the observations one would natu- 
rally expect the galactic disc dynamo to generate this dipole solution. 



26 



which both avoids the flux expulsion problem and has a lower critical 
dynamo number. The fact that the Galaxy possess the more difficult 
quadrupole solution with its larger critical dynamo number seems to 
imply that the field did not originate from a seed field with random 
symmetry, but had an already established quadrupole symmetry and 
a respectable field strength from the beginning. We discuss the flux 
expulsion problem in §9. 

We conclude this section with an estimate of the dynamo number. 
The estimates of the value of ao and /3 given by Parker, and RSS, will 
be presented in the next section. It will be shown that earlier estimates 
converged on a value of a of about 0.1 km/sec. Taking the above val- 
ues of P and h, and Q, 10~^^/sec, we find that D is about unity, 
a value an order of magnitude too small for dynamo amplification of 
the quadrupole field. However, as seen, in §7 the properties inferred 
for the interstellar turbulent motions driving the dynamo are too pes- 
simistic. [See Field (1995)]. During the nineties Ferriere (1992a,b, 
1993a,b) recognized that the observed turbulent motions were mostly 
the intrinsic motions of expanding supernovae remnants and expand- 
ing super bubble shells. These motions have much longer correlation in 
scale and time than that previously assumed for the turbulence, and 
thus lead to larger values for the dynamo number and for the growth 
rates. We discuss the dynamos based on this turbulence in §8. 

7 The magnitudes of a and f3 in the 
interstellar medium 

As mentioned in the previous section, the most popular theory for 
the origin of the Galactic field is based on the mean field dynamo 
theory first developed by Steenbeck, Krause, & Radler (1966). For a 
systematic presentation see Krause and Radler (1980). This theory 
was applied to the Galactic disc dynamo by Vainshtein &: Ruzmaikin 
(1972) in the early seventies. At the same time and completely inde- 
pendently of the work of Steenbeck et al., and of Vainshtein Sz Ruz- 
maikin, Parker developed his own dynamo theory (Parker 1970a, b,c, 
1971a,b,c,d,e, 1973a, 1979), and his own application to the galactic 
dynamo problem. 

The work of Vainshtein and Ruzmaikin is best summarized in the 
book "Magnetic fields in galaxies " written by Ruzmaikin, Shukurov, 
and Sokoloff (1988). which we will refer to in the remainder of this 
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review as RSS. Parker's work is expanded and detailed in his book 
"Cosmical Magnetic Fields ", (1979). The equations that Parker de- 
rived are completely equivalent to those of Steenbeck, Krause and 
Radler, and in his application to the Galactic field, are also equivalent 
to those of RSS. 

The success of the dynamo theory for an origin of the Galactic 
field that is based on the amplification of weak fields during the life of 
the galactic disc, hinges on the values of the coefficients a, (5 and the 
effective Galactic half thickness h. These are encapsulated into the 
dynamo number D, in equation (16). 

Let us assume that when the Galactic disc is first formed, it pos- 
sessed a very weak field with field strength ranging from 10~^^, (which 
can arise by compression from an even weaker field strength 10"^'^) 
generated earlier by a Biermann battery, or by other means. 

Then, to increase this field to its present strength of a few micro- 
gauss and quadrupole symmetry, certain conditions must be satisfied. 
The dynamo number D defined in the previous section must be neg- 
ative and exceed 13 in magnitude by a large enough amount that the 
growth rate is fast enough to increase the Galactic field from its initial 
very small value to its present value of a few microgauss. For exam- 
ple, if the growth time were five hundred million years and the age 
of the Galactic disc ten billion years, then the field must exponenti- 
ate twenty times an increase by a factor of 10^. This would amplify 
an initial field of order 10~^^ gauss to the present observed value. A 
somewhat faster growth than this could amplify the field from 10~^^ 
gauss. The question reduces to whether the interstellar values of the 
coefficients can achieve so fast a growth rate. 

Both Parker (1979) and RSS give rough estimates of the dynamo 
number and when reduced to the same definitions the values based on 
rather diff'erent estimates come out to be remarkably close. 

In this discussion we defer the problem of fiux escape and the 
justification for vacuum boundary conditions to §9. 

The critical dynamo number for any growth at all is D = —13. The 
dimensional growth time, /i^//?, was shown in the last section to be 
five hundred million years. If we take such a growth time as adequate 
for dynamo generation of the field, then we need the dimensionless 
growth rate 7' 1. 

Then from the figure VII. 2 in RSS we see that the required dynamo 
number \s D = —25. If the absolute value of D is smaller, then the 
growth time is longer. 
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Parker estimates D in liis boolc (Parker 1979), making the best 
estimates he can of a and [3. He bases his estimate of a (T in his 
notation) on a model he constructs for a rising convection cell in the 
Galactic disc. He arrives at the estimate 

a K -irem OMLU, (17) 
8 

where for his model e 0.1, and L is the radius or spacing of convec- 
tion cells. For this he takes L lOOpc. This leads to 

a « O.lkm/sec . (18) 

He estimates /3 (which he calls rjr) from another model for diffu- 
sion. He finds that 

P ^ 0.2vL (19) 

where v is the observed velocity of the interstellar clouds which he 
assumes is typical of the observed interstellar turbulent motions. 
Thus, his value for D (which he calls {Kh)^ ) is 

{Khf=D^-{^)ih/LfinTf (20) 

where r = L/v is his estimate of the correlation time. He takes h = 
400ps . His factor of 3/2 is based on assuming Kcplcrian rotation 
of the Galactic disc, and should actually be replaced by unity. This 
equation (after this replacement) leads to his estimate of the dynamo 
number D ~ 4. 

This estimate is based on very rough assumptions, so the fact 
that it is less than the critical number for excitation of the a — fl 
dynamo in the the quadrupole mode is not of too great a concern. It 
is interesting that in spite of the very large and small numbers going 
into the definition of the dynamo number it turns out to be of the 
correct order of magnitude for exciting the dipole mode. 

RSS also made an estimate of the dynamo number. Their estimate 
is based on a less careful estimate of the a — $7 dynamo coefBcients, 
but leads to a value rather close to Parker's estimate. Their estimate 
for a given on page 181 of their book is a = /Q,/h « Ikm/s , but 
apparently they left out the factor of one third in the definition of 
a. Taking the same number L = lOOpc, h = 400pc we see that their 
estimate for a is larger than Parker's estimate of P = .04Ln by a 
factor of about six. (With the factor of one third th is would be twice 
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Parker's value.) Their estimate for /3 is /3 « -^Lv larger than Parker's 
estimate by a factor of five thirds. These values lead to their estimate 
for D m 10. If we restore the missing factor of 1/3 in RSS's estimate 
for a, their estimate reduces to D ^ 3, not far from Parker's estimate 

Both of the two estimates for D are substantially smaller than that 
necessary to produce the present value of the Galactic field starting 
from a reasonable seed value in the lifetime of the Galaxy. 

We have given these estimates in some detail because for a while 
these were the principal estimates for the action of the a — fi Galactic 
dynamo and the astrophysical basis of these estimates can be easily 
understood. If taken seriously, they result in a rather discouraging 
numerical situation for the Galactic dynamo. Clearly, to make the 
dynamo work some way to increase the effective value of a is needed. 

It is interesting that for closed boundary conditions (no flux es- 
cape) the critical dynamo number is about -4 for the dipole mode. 
These boundary conditions correspond to no net flux in the local re- 
gion of the Sun. Since the observations seem to require the nonzero 
flux quadrupole mode, one has a feeling that if the Galactic dynamo 
is the origin of its field, the Galaxy would choose the more difficult 
quadrupole mode for its magnetic field structure. 

It is intriguing and provocative that more recent Faraday rotation 
observations seem to indicate that there actually is some asymmetry 
about the Galactic midplane (Han, 2001, 2002, 2006, Han Manchester, 
Berkhuijsen and Beck 1997) These results are based on Faraday rota- 
tion measurements of pulsars at Galactic latitudes 6 > 10 deg in the 
second and third quadrants. The other measurements still are quite 
symmetric. 

As mentioned in the last section, the remaining option is, of course, 
that the magnetic field is pregalactic, and initially had the quadrupole 
symmetry. 

This was the dynamo situation in 1991. But the whole picture 
was changed substantially by a careful analytic series of excellent pa- 
pers published by Katia Fcrriere (Ferriere 1992, 1993a, b, 1996, 1998, 
Ferriere &: Schmidt 2000). These papers turned the whole situation 
around leading to a much more definitive analysis of the astrophysical 
situation, a more precise evaluation of the dynamo coefficients, and a 
much clearer picture of the helicity of interstellar turbulence. It turns 
out that the turbulence of interest is not that resulting from stirring 
of the interstellar medium by stellar winds and supernova. The mo- 
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tions of importance are those actually present in the supernovae and 
sTipcrbubblc explosions themselves. We discuss Ferriere's analysis in 
the next section. 



8 Ferriere's dynamo theory based on 
supernova and superbubble explosions 



The estimates of dynamo action in the Galactic disc by Parker and 
RSS were derived from what little was known in the seventies and 
eighties about interstellar turbulence, as well as an educated guess as 
to the effect galactic rotation would have in producing kinetic helicity. 
It has always been correctly assumed that the main source of the 
observed turbulent motions is due to supernova explosions and stellar 
winds. However, it was earlier supposed that this turbulence arose 
due to the injection of supernova energy into the interstellar medium, 
stirring it up into turbulent motions. This did not appear to make 
the turbulence strong enough to amplify the Galactic field. 

Ferriere took a more direct approach to the dynamo problem, us- 
ing the actual velocities of the expanding supernova shells during the 
explosions themselves as the mechanism for producing a and /?. That 
is, she asserted that the magnetic field is produced during the actual 
explosions, and much of the random motion we sec in the interstel- 
lar clouds and interpret as fluid turbulence, is the actual expansion 
velocities of the explosions themselves. 

To analyze this, she made use of a blast wave model for supernovae 
in which the interstellar magnetic field and the ambient interstellar 
medium at a particular time are swept up into a thin rapidly expanding 
sheet. She employs a definite expression for r, the radius of this sheet, 
as a function of time 



(Ferriere 1992a, Weaver et al, 1977). She then derives the behavior 

of the plasma motions and the magnetic field evolution in this sheet, 
under the influence of Galactic rotation. As the shell expands, its 
moment of inertia increases by a large factor and the shell tries to come 
to rest in the laboratory frame, thus appearing to rotate backwards in 
the frame of the rotating Galaxy. Since the initial field and the matter 
in which it is embedded were initially rotating with the Galaxy, the 




(21) 
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effect of the Coriolis forces produce helical motions that twist the 
magnetic field out of its original plane. If the field is initially toroidal, 
then this twist partially rotates it into the radial, or poloidal, direction 
and produces the a effect needed to complete the dynamo. 

The actual angle through which the field is twisted is roughly the 
angle that the Galaxy rotates at the radius of the supernova, i.e. the 
angle it would rotate during the lifetime of the supernova, a few hun- 
dred thousand years. This twist is actually quite small, so that a 
perturbation calculation is adequate for the twist of the field. How- 
ever, under the influence of the many supernovae that a given piece 
of the interstellar medium suffers, the twist accumulates. Since it is 
always in the same direction, it accumulates to a finite value in few 
hundred million years. 

Ferriere calculates this effect with considerable precision taking 
into account that, relative to a fixed point, P, the many different 
supernova remnants which overlap this point occur with different di- 
rections and distances from P. Some supernovae occur quite close 
to P, and produce a stronger effect than those that explode farther 
away. It is known that the frequency of supernovae at any given place 
is roughly comparable to the reciprocal lifetime of a typical supernova 
remnant, so that the angular rate of the accumulated twisting of the 
piece of the field line at point P is a finite fraction of the rotation 
of the Galaxy. Thus, it is no longer necessary to make a very rough 
guess as to the value of a. By using her technique the full a tensor 
can be calculated with some precision. Its value depends mainly on 
the frequency and spatial distribution of the supernovae. These are 
reasonably well known, as are the details of the expansion. 

Further, it was possible to calculate the other parameter of the 
dynamo theory, /3. This is obtained by evaluating the mixing of the 
plasma and magnetic field while they are in the expanding shell. 

There is another important parameter she introduces into the dy- 
namo theory, the escape velocity, T4sc of the field lines and plasma. 
Because supernovae occur closer to the midplane of the Galaxy than 
the bulk of the interstellar medium and its magnetic field, the result 
of any explosion on the average is to expel the plasma and field lines 
away from the midplane. Ferriere points out that this average dis- 
placement per unit time is equivalent to Vesci but she keeps it in the 
turbulent term in the dynamo equation. 

A difficulty with her calculation is that she only includes the effect 
of the supernova remnant on the plasma and field when the remnant 
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is expanding, and does not consider the effect on the field when its 
plasma falls back to its original position, after the supernova expansion 
concludes. Her assumption is that the magnetic field would not be 
acted on coherently during the fall back, so that the twisting that 
occurs during the explosion would persist. This again involves the 
physical process of the magnetic flux separating from the plasma. This 
remains a serious problem in her calculations which we also discuss in 
S9. 



A pictorial view of how the supernova explosion actually ampli- 
fies the magnetic field without violating flux conservation is shown 
in Figure 1, (Kulsrud, 2000) The various panels are alternate views 
from the side towards the Galactic center and from above. In panel 
(a), the supernova is shown just before it goes off and an undisturbed 
magnetic line in the toroidal direction is also shown. After the su- 
pernova remnant has expanded to a radius p and the radiative phase 
of the explosion is reached, nearly all of the swept up matter is con- 
tained in a thin shell at its boundary and the magnetic field lines are 
embedded in this shell. This is shown in panel (b). Outside of the 
bubble the field and plasma are undisturbed. Because the matter has 
been shifted from the place it occupied before the explosion where it 
uniformly filled the sphere to the shell, the matter in the shell has 
increased its moment of inertia by a large factor. By conservation 
of angular momentum the supernova bubble decreases its rotational 
velocity in the laboratory frame and rotates backward in the Galactic 
frame, as shown from a top view in panel (c). The piece of the orig- 
inal line in its original plane is shown by the solid line, and the part 
overlapping the supernova bubble is shown by the dotted line. 

The rotation has displaced the footpoints of the line in the radial 
direction in this plane. But this displacement plus the differential 
rotation of the Galaxy displaces the lines in the toroidal direction, 
as shown in another top view in panel (d). Thus, in the original 
horizontal plane the field line has doubled over itself, increasing its 
field strength and the magnetic flux in this plane. On the other hand, 
the field overlapping the supernova has developed negative toroidal 
flux compensating for this increase and still conserving total toroidal 
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Figure 1: The operation of the a — Q dynamo. In (a) and (b) a supernova 
blows the magnetic loop. It is twisted into the poloidal plane by Coriolis 
forces as is seen in the top view of (c). The dashed line represents the upper 
part of the loop. In (d) the lower part of the line is stretched by differential 
rotation. Then the upper part of the line is removed to infinity as in (e) and 
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flux. If the supernova continues to expand indefinitely, or at least the 
flux line continues to rise indefinitely, as shown by the side views in 
panels (e), and (f), the part of the field line on top of the supernova 
bubble would be removed from the Galactic disc and the toroidal flux 
in the disc will have increased. Since negative flux is removed far 
away, it is no longer counted when one describes the Galactic field. 
No Faraday rotation measure passes through this removed flux and 
the rotation measure in the disc is now doubled. In this manner the 
problem of flux conservation is solved. [Note that if the two vertical 
parts of the line in panel (f) should magnetically reconnect then the 
lower part of the reconnected line would return to the disc and cancel 
any net gain in magnetic flux there.] 

The last of the panels in this picture represents how the vacuum 
boundary conditions, assumed in the conventional a — Q dynamo, are 
imagined to act. One can see the difficulty with which this dynamo is 
faced. The question is reduced to whether the flux in the expanding 
plasmas can actually escape from the Galaxy. Generally, except for 
supernovae that occur in the halo, the supernova bubble is slowed 
down below the Galactic escape velocity by gravity and ram pressure 
on the interstellar medium. It stops long before it reaches more than 
a hundred parsecs or so above the plane. Then the danger to the 
dynamo field amplification is that the matter in the supernova shell 
falls back down and any unwanted flux still embedded in it cancels 
out the amplified field. The result in this case is that there is little or 
no amplification at all. 

The separation of flux from plasmas during the expansion is crucial 
to the amount of amplification. If the field is very weak, there is no 
way for the plasma to know it has flux embedded in it, and there will 
be no separation at all. 



To calculate the dynamo action of a supernova Ferriere selects 
some point P at which she calculates the dynamo action. (See Figure 
2.) She then calculates the electric fleld E at this point during the 
time in which the shell passes over it. This is given by 

^ <^v X B 

E = (22) 



35 




Galactic Midplane 



Figure 2: The supernova shell passing point P 

at this point while the thin supernova shell passes over this point. 
(5\- is the expansion velocity which though large she denotes by a 5 
since she regards it as random when summed over many supernovae.) 
She treats the Coriolis forces and the magnetic field twisting as first 
order in a perturbation expansion. To first order she takes B from 
flux freezing. That is she calculates the field as though it were simply 
swept up by the snowplow effect of the supernova shock using the blast 
solution for the shock motion given in equation (21). 

With this model the electric field at the point P is nonzero only 
when the shell passes over it. There is no velocity before the shell 
reaches P, and after it passes, the magnetic field is reduced to zero. 
The time integration of E at P is just such as to produce the field which 
cancels the original field leaving B zero as one expects. Although the 
shell may be very thin, and the time over which the E is nonzero 
very short, the magnetic field in the shell is very large and the actual 
magnitude of E multiplied by the time when it is nonzero is finite and 
independent of this thinness. 

To properly include the twisted part of the changed magnetic field 
Ferriere evaluates the first order velocity 5v in the shell, taking into 
account the effect of the Coriolis force on the zero order velocity, and 
5B, the first order change of the magnetic field produced by the (5v. 

Substituting her result for Sv and (5B into the expression for E 
given by equation (22), integrating it over time at the fixed P, sum- 
ming this result over all supernovae that interact with the point P 
during some long time t, and then dividing by t, she finds the long 
time average of E at P. She assumes supernovae occur randomly with 
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a distribution that depends only on z, the height of above the plane. 
The resulting value she finds for E gives her the mean value for a 
and P to use in the dynamo equation to calculate the evolution of the 
mean field B. See equation (11) of §6. 

The result is a tensor for a which depends on the original orienta- 
tion of B before she started her averaging over supernovae. If B is in 
the toroidal direction, i.e. the same direction as the Galactic rotation, 
then a is a scalar, the same scalar originally considered by Parker, 
RSS, and others. However, in general, B will have other components 
and the tensor character of a becomes important. (Ferriere, 1998). 

In addition to a she also obtains expressions for the /? tensor and 
for Vgsc; which is a vector in the z direction. 

Using the same technique, Ferriere also considered the contribution 
of superbubblcs to the Galactic value of a. A superbubble is the 
result of many supernovae going off effectively simultaneously. The 
simultaneous explosion of many supernovae at the same time and place 
arises from the formation of stars close together in clusters at nearly 
the same time. The more massive stars exhaust their nuclear fuel at 
nearly the same time and one finds that, during the lifetime of a single 
supernova, many others explode. The gaseous remnants overlap and 
form a low density bubble of high pressure gas. This drives a shock 
that arises from the energy input of all the supernovae, rather than 
from the kinetic energy of the individual explosions. 

As far as driving the dynamo effects the superbubble acts just as a 
single supernova does. In the thin shell, driven outward at high veloc- 
ity, the Coriolis forces act to twist the lines. However, a superbubble is 
much bigger than a single supernova shock and lasts much longer. As 
a result, the effect on the dynamo is much greater than that of individ- 
ual supernovae, even when one takes account of the smaller number 
of superbubblcs than supernovae. The radial time dependence of the 
superbubble shock which Ferriere employs is (Ferriere 1992a, McCray 
k Kafatos, 1987) 



where L38 is the luminosity in units of 10^® ergs/s and is the cylin- 
drical radius. 

Comparing the times and distances in this equation with those 
for the single supernova, one can sec that the superbubblcs extend 
to much greater distances and last a great deal longer. The larger 
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distances and times allow the galactic Coriolis forces to act longer, 
and produce a larger effect in the superbubble shells than in single 
supernova shells. This is compensated somewhat by the fact that 
there are typically thirty times more supernovae than superbubbles. 
However, even taking this into account Ferriere found that the effect 
of all the superbubbles was larger than that of all the supernovae by 
a factor of order seven, even though the total number of supernovae 
in superbubbles is comparable with their total number of supernova. 

Making use of these calculations, Ferriere developed the Galactic 
dynamo to a new level of precision. She established several new results 
in these papers. One result was that the a, in the alpha effect is a 
tensor unless the initial Galactic magnetic field is purely in the toroidal 
direction. She calculated the P tensor using the same procedures. In 
addition, she found a new systematic velocity Vesc that represents the 
rate at which the magnetic field lines rise above the midplane of the 
Galaxy (during the expansion phase of the bubbles). Ignoring any 
fallback of the magnetic flux, she took this velocity to represent the 
rate of flux loss to the disc. 

There is one correction that should be made to her calculations. 
She calculates the mean a effect at a fixed point. That is to say, 
she assumes that the a tensor which she calculates gives the rate of 
evolution of the poloidal field at this fixed point. But there is a term 
missing in her calculation that when introduced implies that her value 
of a actually gives the evolution following the expansion of the shell. 
This seems intuitively clear when one notes that, after the shell passes 
point P, there arc no field lines left at P. All the field lines have been 
assumed to be swept up by the shell. This extra term disappears if one 
carries out the calculation in the moving frame and expresses a as a 
function of the Lagrangian coordinate, the fluid element. Her method 
of calculation working with the averaged time integral of the electric 
field, which she calls the E.M.F., tends to hide this extra term. 

Correcting for this by making this simple transformation from an 
Eulerian independent variable to the Lagrangian variable should make 
her calculation more correct. 

It is interesting that she calculates the actual change of the field in 
Lagrangian coordinates, and then uses this changed field to calculate 
the dynamo coefficients a and /? which are subsequently used to calcu- 
late the mean change in B. This roundabout procedure is an attempt 
to relate her calculation to the standard a — Q mean field dynamo 
theory. However, the approach from supernova and superbubble to 



38 



the dynamo does not quite fit the assumptions of the dynamo theory 
which assume that the turbulence scales are small compared to the 
large scales on which the mean field is to be calculated. The super- 
nova and superbubble remnants are on a scale of the total thickness 
of the Galactic disk. It would be neater if she just stopped with the 
changed fields in any given Lagrangian element of plasma and directly 
averaged the changes occurring in a given plasma element rather than 
referring them to the dynamo coefficients. 



9 The validity of the vacuum bound- 
ary conditions 



As emphasized in §6 and §7, boundary conditions at z = zth, are 
needed to compute the rate of growth of the dynamo mode. To see 
that these conditions are directly related to the change of flux in the 
disc for the quadrupole eigenfunction, integrate the first equation of 
(13) from -h to +h. This yields 



where Bj., Bq are assumed symmetric. We see that the first term 
—/3dBr/dz is the rate of escape of the Br flux through the boundary 
while from the zero boundary conditions on Br and Bg, the second 
term vanishes. Although the flux lines of the mean field are not frozen 
in the plasma its global fiux inside the disc still satisfies the conserva- 
tion relation (24). The Bq flux is not per se conserved, since it grows 
as the Br flux is wound up, but this does not increase the number of 
lines of force. 

As discussed in §6, the boundary conditions customarily assumed 
at ±/t are that the fields vanish, and this is the natural condition if the 
flux which reaches the boundary, escapes instantaneously. Physically 
this is difficult to imagine because, flux freezing implies that any flux 
that escapes is embedded in interstellar matter. Thus, for the dynamo 
to work, a large portion of the interstellar medium must be removed 
at every efolding of the magnetic field. [However, see Brandenburg, 
Moss & Shukurov, (1995), Doblcr, Poezd, & Shukurov, (1996) and 
Shukurov, (2004), where this problem is discussed.] Let a fraction 
/ of the entire interstellar medium be removed for every exponential 
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increase of the magnetic field. Then, if the field strength increases 
from Bq to Bi, the mass of the interstellar medium decreases from 
Mo to Ml = Mo{Bi/Bo)-f. To raise the field from Bq = IQ-^^ 
gauss to Bi = 10~^ gauss with / = 1/3 would require Mq ~ 2000Mi, 
which seems untenable, and inconsistent with the progressive chemical 
enrichment of the interstellar medium over time observed from stellar 
composition measurements. 

These numbers give one pause when considering these assumed 
boundary conditions. 

Further, the energy to remove this amount of mass to infinity re- 
quires l/2Ug^^ per unit mass, where v^sc ~ 500km/s is the escape 
velocity (Smith et al. 2007). This gives of order 10^^ ergs /g. Taking 
Ml ~ 10^ solar masses we find an energy requirement of about 10^° 
ergs. Thus, an amount of order the total luminosity emitted over the 
disc life, is needed to amplify the magnetic field by the galactic dy- 
namo under the above assumptions. Clearly, a better way is needed 
to remove flux. 

A more analytic treatment of the boundary conditions is discussed 
in the review article of Beck et al 1996. They state that the turbulent 
resistivity /3 is much larger by a factor of say in the halo, (the region 
outside the disc \z\ > h) larger than in the disc, and that a vanishes 
in the halo. In this case, the components of the halo magnetic field, 
in the dimensionless units of equation (14) of §6, are proportional to 
e~-^^ and the boundary conditions at 2;' = ±1 should be that the field 
there is smaller than the average field by about i^l/A). For large 
A the boundary conditions become essentially the vacuum condition 
B = 0. 

To justify this large increase in the halo value of /3, Beck et al. 
quote results from Poezd, Shukurov k. Sokoloff (1993). These authors 

estimate (3 from measurements of the random velocities in the halo 
made by Kulkarni &; Fich (1985). The results were %) 10 km/sec, 
and the correlation length £ fs 100 pc. This increases /3 by a factor of 
one hundred above its disc value. But these measurements are made 
on the bare halo at a time when there is no escaping flux or mass. 
Surely if there was a large amount of mass diffusing outwards, the 
velocity fluctuations would be much smaller due to the large amount 
of energy involved in lifting the mass-loaded flux outward. 

Beck et. al. (1996) also appeal to arguments given by Brandenburg 
et al (1995), who attempt to treat the actual outward escape of flux by 
invoking reconnection of it with the halo field. Again, the reconnection 
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for this merging is based on replacing the normal resistivity 77 by the 
very large turbulent resistivity (3. 

Now it is not clear that there is enough magnetic flux in the halo 
field to match the expelled flux and accommodate this reconnection. 
Further, /? is not a true resistivity and can only mix the field lines of 
different sign. The field lines rising up from the disc would still be 
loaded with the heavy interstellar mass while the other lines would be 
much lighter. Then even after the lines axe mixed, the gravitational 
field could separate them allowing the original field lines to sink back 
and cancel any gain from the dynamo. 

It must be borne in mind that, when trying to arrive at the origin 
of the Galactic magnetic field, we are primarily interested in building 
up the field from an extremely weak initial value. When the field is 
weak, such mixing is unlikely to occur. Also any true reconnection 
due to real resistivity would be too slow to be of interest. 

A third possibility to justify the boundary conditions, at least 
physically, is to imagine that the escaping flux is lifted in huge arcs 
such as occur in superbubbles, and to have the mass slide down the 
flux tubes. This would leave the field lines at the tops of the arches 
unloaded and lighter so that this piece of them might escape by buoy- 
ancy or cosmic ray pressure (Parker 1992, Hanasz et al 2004). How- 
ever, again this can not work in the early dynamo stages when the 
field is weak because the magnetic field can have no dynamic conse- 
quence; the fluid flows in a direction independent of the magnetic field 
and probably just falls straight down. 

It also docs not seem to work when the field is stronger, at least 
for the superbubble case, because the mass at the top is on a nearly 
horizontal line, and the rate of flow down the field lines is too slow 
to release much flux. This is demonstrated in the paper of Rafikov 
& Kulsrud (2000). Further, the presence of cosmic rays on the field 
lines in this latter case also inhibits the downflow along the lines. 
This is because an Alfven wave instability is produced by any relative 
motion of the cosmic rays and the matter exceeding the Alfven speed. 
(Rafikov & Kulsrud 2000, Kulsrud & Pcarcc 1969.) This instability 
couples the parallel motion of the plasma and the cosmic rays and has 
the result that the cosmic ray pressure stops the downsliding. 

The most reasonable way for the flux to escape and properly com- 
plete the ail dynamo theory would be to have the superbubbles blow 
the mass entirely out of the Galaxy. But, in our Galaxy, except for 
the rare case where a superbubble occurs high in the halo, the su- 
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perbubble shell is slowed down by the ram pressure of the swept up 
material. Then the gravitational downpuU is enough to stop it before 
it has gone far enough to escape from the Galaxy. This is found to 
be the case in the many numerical simulations (MacLow & McCray 
1988). 

However, as emphasized by Heiles (1990) the number of supernovae 
in a superbubble assumed in these simulations is taken as the average 
number (thirty) and further the superbubble is taken to start at the 
midplane. Superbubbles which contain more than the average number 
of supernova and which start at a position that is high enough above 
the central plane of the disk that lies above the bulk of the interstellar 
medium will lead to a blowout that expels their mass entirely out of 
the galaxy, (See in figure If.) These events which are rare, should 
occur independently of the magnetic field strength and happen at any 
time during the life of the disc, even when the magnetic field is weak. 

As long as the vertical field lines do not reconnect some net flux 
expulsion will occur and make the galactic dynamo viable although 
if these events are rare the dynamo will be weak. One cannot ig- 
nore this possibility. However, if this is the case, the dynamo would 
be extremely inhomogeneous. The resulting field lines would be dis- 
continuous as viewed inside the disc (but actually connected through 
intergalactic space). Such a configuration has not been considered 
in interpretations of the galactic field. Whether it could be consistent 
with observations of the interstellar field is unclear. However, it would 
certainly change the general picture of the galactic field. 

For smaller galaxies in which the gravitational field is weaker or 
in starburst galaxies where the superbubbles are stronger the flux 

expulsion problem is eased and the a O dynamo actually could 

more easily produce magnetic fields. 

The OL 17 dynamo is a beautiful mechanism. It only appears 

to be limited by these difficult boundary conditions issues. However, 
they arev very hard to get around when the magnetic field is very 
weak. In fact, it does work quite well for the dipole symmetry where 
the total flux is zero and does not need to change. If the q — dynamo 
is the main agent for growing the magnetic field strength, it seems very 
surprising that the Galactic symmetry is not odd. 

It is entirely possible that a reasonable way to separate out the un- 
wanted flux may emerge and resolve the Galactic dynamo problems. 
This would go a long way toward establishing an origin for cosmic 
magnetic fields. Any indirect evidence as to the state of the Galactic 
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field at the inception of the disc would also go far toward resolving 
whether the a Q dynamo is the primary source of the field. In- 
deed, there is actually some evidence which, while still uncertain, does 
bear on the problem. This evidence, discussed in §2, comes from the 
abundance of the light elements Li, Be, and B in very old stars. 

10 Arguments against a primordial ori- 
gin 

We have raised some difficulties with the galactic disc dynamo as the 
sole origin of the Galactic magnetic field. There are also a number of 
arguments that have been raised against a pregalactic origin. These 
have been expressed by Woltjer (1969), Parker (1973b, 1979 p. 519 
and p. 525), and Rosnesr & Ducca, (1988). [But see Kulsrud (1990), 
Howard and Kulsrud (1997) and Kulsrud (2006).] 

The mechanisms for the production of a pregalactic field have been 
developed with much less precision than those for the galactic dynamo. 
The plasma physics of the very early universe is only partially known. 
In addition, events that occur very early on have a coherence scale 
smaller than the Hubble radius at the corresponding redshift z. In 
fact, the corresponding comoving scale (the present scale of the cosmos 
that had the Hubble radius at red shift z) is equal to the present 
Hubble radius times (1 + z)"^/^. The observed scale of the Galactic 
field is at least 1 kpc which corresponds to a region of the universe of 
about 10 kpc before collapse to the Galaxy. Thus, any magnetic field 
generated at a red shift z > 10^^, corresponding to a temperature 10^ 
electron volts, would lead to a field whose present scale is smaller than 
1 kpc in the Galaxy. After the time corresponding to this temperature, 
which is the electron - neutrino decoupling temperature, the plasma 
physics seems well understood and there are few grounds for believing 
that anything mysterious could generate a magnetic field. 

A very weak field with strength of order 10~^° gauss could rea- 
sonably be generated by action of the Biermann battery mechanism 
discussed in the next section, (Biermann 1950), or by interaction of 
the electrons with the CMB photons (Mishustin & Ruzmaikin 1972). 
But such a field could only serve as a seed field, and this gis generally 
accepted as a reasonable seed field for galactic dynamos. 

One possibility that has been taken seriously for a pregalactic ori- 
gin, is a field generated at moderate redshift by turbulence that occurs 
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during structure formation, say in the protogalaxy, (Pudritz &; Silk, 
1989, Kulsrud, Ccn, Ostrikcr & Ryu, 1997). This possibility will be 
discussed in §12. Another possibility that is seriously considered is 
that the magnetic field is somehow generated when a large radio jet 
is formed, say by dynamo generation by a black hole at its center. 
Magnetic fields of substantial strength are presently observed in such 
giant radio jets. There are even some Faraday rotation measurements 
of such fields that indicate the magnetic flux is large enough that if 
spread out throughout the universe it would provide the total fiux 
currently observed in all the galaxies. If this should happen before 
galaxies form, then the field is already present in the plasma out of 
which the galaxies form, and would thus have a pregalactic origin 
(Rees & Seti, 1968, Daly & Loeb, 1990, Furlanetto & Loeb, 2001). 
So far, the details of this filling of the universe with fiux have not 
been worked out with any precision, but the presence of this amount 
of flux, if supported by the Faraday rotation measurements, is very 
suggestive as an origin for pregalactic fields, (Colgate &; Li,2000, Li at 
al. 2006a,b). 

Before taking up the discussion of the first of these proposals for 
a pregalactic origin let us discuss the arguments that have been made 
against any primordial origin as opposed to a galactic dynamo origin, 
(see Kulsrud, 1990.) 

There are essentially two arguments against a pregalactic magnetic 
field: (1) the winding up argument and (2) the escape of any magnetic 
field from the Galactic disc in a short time by turbulent diffusion. 

One can explain the winding up argument as follows: Suppose one 
starts out with a uniform unidirectional field Bq in the x direction that 
extends over the entire Galactic disc. Then in cylindrical coordinate 
the radial component of the field is Br = Bq cos 6. Now, at each radial 
circle r of the Galaxy this radial field is rotated at a rate ^{r), so after 
a time t the field is rotated by an angle Qt, and 

Br{r,9,t) = Bocos{9-n{r)t). (25) 

By flux freezing, this field develops a toroidal component 

Be = BJr,e,t)xr^t. (26) 
dr 

This is a slowly varying function of r times Br- But, if we look at 

the variation of Br with r at fixed 9 and t, we see that Br oscillates 
in r with a wavelength equal to Ar = 2Tr /{t\di},/dr\). Now, near the 
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sun, J7r is a constant, so d^/dr = —Q/r and Ar/r = 2Tr/flt «i 1/50. 
Thus, the toroidal field would change sign every 100 pc. Parker states 
that this is contrary to observations, and clearly such a field would 
produce a zero average Faraday rotation measure. 

This winding up argument is one of the basic arguments against a 
primordial field origin. However, we are not sure that the correlation 
distance of the Galactic field is 1 kpc. We can only be sure that if 
it does reverse rapidly, it does not average out over such distances, 
since there are rotation measurements of pulsars much further away 
than this, (Rand & Kulkarni 1989). (The one-kiloparsec correlation 
distance is derived from models that take into account the possible 
reversal of the field between the spiral arms that have this spacing, 
but ignore the models with short reversals.) 

On the other hand, the resolution of pulsar measurements is not as 
fine as 100 pc so that it is entirely possible that the field does reverse 
over 100 pc but does not average out. In fact, if the initial Bq field 
were not uniform, but varied by a factor of two over the diameter of 
the Galaxy, then an inspection of the above argument shows that the 
toroidal component of the field Bg would still reverse with a scale of 
100, but the toroidal field would one sign would be twice as large as 
that of the other sign and the Faraday rotation measurements would 
not average out. (See Howard & Kulsrud, 1997 and §15 for a similar 
result that would occur with ambipolar diffusion). It would be an in- 
teresting challenge to find out whether this is the case, or whether the 
field is actually of constant sign over such small distances. The field 
of our example could account for the puzzling fact that the fluctua- 
tions in synchrotron polarization direction are of different amplitude 
than the fluctuations associated with Faraday rotation measurements 
(Zweibel & Heiles 1997, Brown & Taylor 2001). This discussion is 
in no way able to establish that the Galactic flcld docs indeed reverse 
over a small length. But the possibility of it does weaken the argument 
against a primordial origin. 

The second argument, the escape argument, described by Parker 
on page 522 of his book, concerns the turbulent resistive decay of a 
pregalactic magnetic field. The decay rate by ordinary resistivity is 
very slow and can be disregarded. However, if one replaces the ordi- 
nary resistivity by the turbulent resistivity, f3 lO^^cm^ /sec, then 
the decay time is 300 million years and Parker argues that a primordial 
field would not survive for the lifetime of the Galaxy. However, we see 
in §9 that although this diffusion can mix the fields in the disc, they 
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cannot remove it from the disc unless the field is strong enough that 
it is able to separate itself physically from the interstellar material. 
(Without the separation, the diffusion mechanism would have to lift 
the interstellar medium entirely out of the Galaxy which we showed 
was energetically very difficult.) 

Thus, neither of these arguments can definitively rule out the possi- 
bility that there was an initial large scale magnetic field of pregalactic 
origin. 

In §15 we describe a modification of such a model in which ambipo- 
lar diffusion plays an important role and an initial field stronger than 
10~^ gauss could be brought up to full strength without any dynamo 
action at all. However, total neglect of the a — Q dynamo is, of course, 
not realistic (Kulsrud 1990, Howard & Kulsrud 1997, Kulsrud 2006). 

There is, an additional consideration that depends more on our 
current understanding of galaxy evolution. Our Galaxy currently has 
of order 10^ Mq of interstellar material. There axe many sources 
and sinks of gas: stellar winds, planetary nebula ejection, supernova 
explosions, infall of intergalactic clouds and small galaxies, star forma- 
tion, and a possible galactic wind. Estimating the rates for all these 
processes leads to the conclusion that the interstellar gas is replaced 
mainly through cycling in and out of stars, on 10^ year timescales, 
about 1/10 the age of the oldest stars in the Galaxy. If the magnetic 
field is in a steady state, dynamo processes must operate over that 
time to act on the new material. Of course, the new material may 
already by magnetized at a level far above the estimated pregalactic 
field strength, although it is unlikely to have the correct topology. An- 
other alternative is that the already existing magnetic field somehow 
mixes into the new material, possibly by ambipolar diffusion (Heitsch 
et al 2004). 

11 Seed fields 

It is reasonable to assume that there was a time before which there 
were absolutely no magnetic fields in the universe. If flux freezing were 
an exact constraint, there would be no fields after this time either. 
This also follows directly from the magnetic differential equation (4) 
of §4. How do magnetic fields get started in the first place? 

The most popular mechanism is the Biermann battery (Biermann, 
1950). A magnetic field can arise from extra terms in Ohm's law which 
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are not included in equation (1). 
These terms arise as follows: 

First, we rederive the magnetic differential equation keeping the 
additional terms in Ohm's law. In fact, Ohm's law can be considered 
to be simply the equation of motion of the electron fluid (Spitzer 1962), 

nem— = -ngC (E + Vg x B) - Vpe + neUig + Fgi, (27) 

where Fgi is the electron-ion frictional force, which is related to the r/j 
term. We drop this term. We also drop the inertial term nem{dvc/dt) 
and the gravitational term ngW-g because of the smallness of the elec- 
tron mass density. Dividing the resulting equation by UgC we get 

C UeC 

To see physically how the Biermann battery works it must be real- 
ized that (in the absence of a magnetic field) any accumulation of the 
electron pressure at any point naturally leads to a tendency for the 
electrons to leave this point. But this automatically leads to a charge 
imbalance producing an electrostatic field, and this field is just strong 
enough to resist this tendency. (The charge imbalance to produce such 
a field is generally minute.) This field is 

E = (29) 
-rieC 

E depends on pg and Ug the electron pressure and density. If Ug 
is a constant in space, the electric field is curl free and purely elec- 
trostatic. Thus, it will drive no current just as a battery with its 
terminals unattached to wires produces no current. However, if Ug is 
not constant, E may not be curl free and there can be a path in the 
plasma around which there would be a potential drop. Thus, a cur- 
rent would flow along this path and produce a changing magnetic field 
which would in turn induce an electric field to balance this potential 
drop. The key to the importance of this term is its dependence on Ug 
in addition to its proportionality to Vpe which enables it to act as a 
battery. 

Now, the resulting current is known to be tiny, so Vg Vj to a 
very high degree of approximation. Further, even if the plasma is 
partially ionized, the electron temperature is kept close to the neutral 
temperature by collisions so that Pe/ng = p/n{l + x) = ^p/ + x); 
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where % is the degree of ionization, which for the moment we take to 
be constant. Thus, 



E + = -— ^. 30 

c e{l + x) P 

Finally, taking the curl of this equation and combining it with 
the induction equation, one gets the modified magnetic differential 
equation, 

— = VxvxB + ^„ ^ — -. 31 

dt ^ ' e(l + x) 

It is well known that the vorticity a; = V x v satisfies a very similar 
equation 

V X (v X w) + vV^u. (32) 



_ = V X (v X w) „ ^ "^2- 

ot p^ 



Thus, as long as viscosity is small, these equations are identical up to 
the factor —e/Mc{l + x)- 

Therefore, eB/Mc = and —{l + x}'^ satisfy the same equations. 
It is reasonable to assume that up until a certain time vorticity is 
also zero. Thus, we have that the magnetic field and the vorticity 
satisfy the same zero initial conditions. (fJ = a; = ). Because the 
differential equations have a unique solution we find that 

fj = ^ (33) 

is valid at least until the viscous dissipation term (or the resistivity 
term) becomes important. 

Thus, when vorticity finally starts to grow, the magnetic field will 
grow with it. However, the resulting magnetic field strength due to this 
term is always tiny. To get an idea of its strength in a gravitational 
forming structure assume that its rotational energy is of the order 
of its gravitational energy, and thus the rate of rotation should be 
comparable to the free fall time, 1/ ^/4TrGp. G is the gravitational 
constant and p the density. Now, e/Mc = 10^ sec/gauss, so we have 

B « 10-V47rGp « 10-^^^/^. (34) 

where n is the hydrogen number density cm~^. Thus, for typical 
cosmic densities, the resulting magnetic field is well below 10~^^ gauss. 

Actually, equation (31) for the evolution of the magnetic field only 
holds as long as the viscosity term plays no role in the evolution of 
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uj; that is, as long as the vorticity itself does not saturate. But when 
the vorticity becomes strong there is Kolmogorov nonlinear coupling 
to small scales, establishing a cascade which is terminated by viscos- 
ity acting on small scales. Effectively, vorticity on all scales reaches 
saturation, even though the effective viscosity on large scales is small. 
Thus, the Biermann battery continues to operate, increasing B above 
the value given equation (33). 

If one examines equation (31), which gives the evolution of B, we 
see that there are two limits. When the magnetic field is extremely 
small, the first term on the right, the dynamo term, is negligible and 
the field grows only by the second term, the Biermann battery term. 
One can estimate its rate of increase by assuming there is a finite 
angle between Vp and Vp, and taking the scale size of both p and p 
to be L. Then the Biermann battery term is of order of magnitude 
J) I pi?' ~ v1/L?^ where Vg is the speed of sound. Thus, increases 
roughly as {vs/L)'^t. 

However, when B becomes larger, the dynamo term becomes im- 
portant, i.e. when vQ/L ~ vl/L?. If, for simplicity, we take v ~ Vg-, 
then this happens when Ug/O ~ L. But Vg/^ is the ion gyration ra- 
dius Tj, so the dynamo term becomes significant when decreases to 
become comparable to L. (Initially, of course, the gyration radius is 
infinite.) After this time the Biermann battery continues to operate, 
assuming that there is still a finite angle between Vp and Vp. But it 
is quickly outstripped by the dynamo term, which increases the mag- 
netic field strength B exponentially. This is because the dynamo term 
is proportional to B. 

Recall that we derived equation (31) assuming the ionization frac- 
tion X is constant. Following the birth of the first massive stars and 
quasars, the Universe was pervaded by a network of ionization fronts 
that spread out from sources of ionizing radiation. During this so- 
called epoch of reionization {z ~ 10; Spergel et al 2007), the low 
density gas was ionized first, leaving behind pockets of denser gas 
that was ionized later. Large scale temperature and density gradients 
in the dense gas interacted with the large electron pressure gradient 
in the front itself, allowing the Biermann battery to operate more 
efficiently. Gnedin et al. (2000) showed that the magnetic fields gen- 
erated in this fashion are somewhat larger than the fields generated in 
vortex sheets (~ 10~^^ — 10~^®G), but are still far too weak to explain 
galactic fields without further substantial amplification. 

The rough assumptions going into these relations cannot be taken 
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too seriously, but they do give a simple picture of what goes on in the 
initial build up of seed fields. To gain more precision concerning the 
dynamo action, we need to specify the velocity field in more detail. 
This is done in the next section, in which the assumption is made 
that the velocity field has a turbulent spectrum. This appears to be 
the case in numerical simulations of the formation of gravitational 
structures, primarily in protogalaxies. (Ryu et al. 1994, Kang et al 
1994, Kulsrud et al, 1997) 

We have devoted considerable space to the Biermann battery be- 
cause it can be treated rather definitively. There are other proposed 
origins for seed fields. The earliest is due to Harrison who invoked 
rotating structures in the early universe (Harrison, 1970). He assumes 
that a blob of plasma rotates in the presence of the cosmic background 
radiation, (and assumes the radiation is non rotating because for suf- 
ficiently large z its mass is larger than the plasma mass.) Then he 
points out that the rotation of the electrons in the blob will be slowed 
down by Thompson scattering off the non rotating radiation field, 
while the ions rotate without slowing down. The resulting current 
will produce a weak magnetic seed field. The seed field will increase 
with time inducing a back electric field which will tend to keep the 
electrons rotating. The balance between these two forces determines 
the actual rate of increase of the field. Under most reasonable assump- 
tions, the resulting field strength is a little smaller, but of the same 
order of magnitude as that of the Biermann battery, < 10~^^ gauss. 

Since the Biermann battery mechanism produces an increasing 
magnetic field, there must also be a back E.M.F. resisting this growth. 
But in the case of the Biermann battery, this back E.M.F. must be 
there to balance the Biermann term in Ohm's law (equation (28)) in 
the first place. One could consider the process in the reverse order, 
for it is clear that the magnetic field must grow to induce an electric 
field to balance this term. 

Another origin for seed fields has been proposed by Martin Rees, 
(Rees, 1987, 1994, 2005, 2006). He suggests that supernova remnants 
similar to the Crab Nebula have fields of order 10~^ gauss due to 
the wind up of the central pulsar. The fields are coherent, but contain 
fiuxes of opposite signs on each side of the remnant. He estimates that 
there could be AT 10^ such remnants in the early stage of the galaxy. 
Of course, the fields from these remnants, when they expand enough 
will overlap and cancel so that the resulting fields will be proportional 
to A^^. But, because they may interact, x may not be 1/2, but could be 
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as small as 1/3. Even so, the resulting fields would be in the range of 
rather substantial seed field! However, one might worry 
that the rms field would be larger than this, perhaps larger by a factor 
of N^~^ unless magnetic reconnection or some other process reduces 
it. The details of this interesting suggestion have not been carried out 
since they depend on plasma processes such as magnetic reconnection 
which are very little understood. One should compare this suggestion 
with the suggestion mentioned in §16, that primordial field can arise 
from radio jets with large magnetic fields of alternate signs and these 
can seed the universe with magnetic fields prior to galactic formation. 
This seeding involves very similar problems of cancellation. 

In any event, there seems to be little difficulty with producing an 
initial field whose strength is of order 10~^° gauss, and that can serve 
as a seed field for dynamo amplification to magnetic fields of larger 
strength. 

12 A protogalactic theory for magnetic 
field generation 

Since there are a number of difficulties with the disc-dynamo origin 
for galactic fields, it is of interest to consider the possibility that the 
field could be generated during the time in which the galaxy itself is 
formed. Such a theory was first proposed by Pudritz and Silk(1989). It 
was later developed in some detail by Kulsrud and coworkers (Kulsrud 
et al, 1997). This latter theory is based on the observation that there 
must be considerable turbulence generated during the initial collapse 
of the cosmic plasma to form the protogalaxy. This turbulence is gen- 
erated by the shocks that result from the steepening of the cosmic fluid 
as the gravitational instabilities become nonlinear. These shocks heat 
the plasma to a Jeans temperature that stops the collapse and leads 
to a temporary virialized state. At the same time the shocks, which 
are generally of finite extent, produce shear velocities that evolve into 
Kolmogorov turbulence. This has been established in numerical sim- 
ulations of structure formation (Ryu et al, 1993, Kang et al 1994). 
The turbulence in these simulations has a Kolmogorov spectrum. Its 
largest eddy is comparable to the size of the protogalaxy, and its total 
energy is comparable to that of the total thermal and gravitational 
energies of the protogalaxy, (Chandrasekhar, 1949). 

Under these conditions a very weak magnetic field will grow be- 
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cause the pressure and density gradients produced by the shocks are 
not parallel. This development was verified by integrating the Bier- 
mann magnetic differential equation (31) of the last section, taking 
the velocities directly from the numerical simulations (Kulsrud et al 
1997). A magnetic field of order 10~^^ gauss develops in this simula- 
tion and is found to satisfy equation (33) at nearly every pixel of the 
numerical grid. However, the expected dynamo amplification is not 
found because this amplification is primarily produced by the smallest 
turbulence eddies whose scale is much smaller than the grid size of the 
simulation. These smaller eddies dominate because their rate of turn 
over in Kolmogorov turbulence is faster than the larger eddies. 

The rate of dynamo growth of the magnetic field is approximately 
their turnover rate. (See the next section.) However, Kulsrud and 
his coworkers used Kolmogorov theory together with the theory of 
magnetic field generation developed by Kulsrud &: Anderson (1992) 
to estimate the rate of growth produced by the smallest eddy. (Kul- 
srud et al. 1997). They find that the growth is much faster than 
the rate of collapse of the protogalaxy In fact, the magnetic en- 
ergy doubles roughly one hundred times in this collapse time. (The 
Kulsrud- Anderson theory is discussed in the next section.) 

During this time the magnetic energy saturates and reaches equipar- 
tition with the energy of the smallest Kolmogorov eddy. Further, the 
Kulsrud — Anderson theory shows that the magnetic field has a scale 
much smaller than that of the smallest eddy. If this were the end of the 
turbulent dynamo process, the resulting field would not be of much 
interest since the size of the smallest eddy is a thousand times smaller 
than the size of the protogalaxy (estimated to be approximately 100 
kpc). After the collapse of the protogalaxy, the field would still be 
coherent on a scale more than a thousand times smaller than the ra- 
dius of the Galactic disk. However, even though the magnetic field 
has such a small scale its energy will be that of the smallest eddy. 
This is of the order of one per cent of the total binding energy of the 
protogalaxy. Such an energy is not negligible, and is strong enough to 
affect the early formation of stars. 

Due to the limited spatial resolution of the simulation, it is impos- 
sible to numerically study the evolution of the field produced by the 
entire turbulent spectrum. In order to form a large scale field by this 
process, there has to be an inverse cascade of magnetic energy to the 
larger scales. Whether this can in fact occur is still an open question 
in dynamo theory. 
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However, an effective inverse cascade can occur by direct amplifi- 
cation of the seed fields by the larger eddies, (see Kulsrud 2005, 2006). 

To predict how the field develops one needs to gain a deeper under- 
standing of the behavior of magnetic fields in the presence of strong 
turbulence, a difficult problem of plasma physics. Schekochihin, Cow- 
ley and their colleagues have published a number of theoretical and 
numerical papers discussing problems concerning the structure and 
saturation of these small scale fields. They consider this problem 
both in the context of protogalaxies, and in the intracluster medium 
in clusters of galaxies. (Schekochihin & Cowley 2004, 2006 Schekochi- 
hin et al 2004, 2005a,b). Their results are discussed in §14. In his 
book, Kulsrud makes an estimate of this behavior (see Kulsrud 2005, 
2006). 

When one considers the small scale turbulent eddies even in a 
rotating body, one finds they have very little helicity. The behavior of 
magnetic fields in non helical turbulence is to be compared with their 
behavior in the a — Q, dynamo theory which is strongly dependent on 
helical motions. The latter theory, combined with proper boundary 
conditions, actually generates net magnetic flux and a mean magnetic 
field in the entire region of interest, while the former theory produces 
fields that are coherent on the scale of the relevant eddy. However, 
the a—i} theory is usually developed in a thin disc and the horizontal 
coherence length of the magnetic field is usually assumed to be a few 
times the disc thickness. 

The radial coherence length of the present Galactic field inferred 
from modeling pulsar rotation measures is at least a kiloparsec, but 
not necessarily larger. Indeed, it seems to reverse in different spiral 
arms. If we compare the distance between the spiral arms with the size 
of turbulent eddies, inferred to exist from the numerical simulations 
of protogalaxies, we see that, if eddies as large as one tenth of the 
protogalaxy generate the main fields, then the resulting fields, after 
collapse, will have a sufficiently large coherence length to provide a 
reasonable origin for the pregalactic field. 

Further, even if the protogalaxy provides a smaller scale field it 
could be large enough for the mean field disc dynamo to work properly 
and develop larger coherence lengths. As emphasized before, the field 
might be large enough to provide the escape of flux needed to validate 
the required boundary conditions. 

This conclusion for definiteness is based on the assumption of top 
down formation of the protogalaxy. It is not drastically modified even 
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if our Galaxy were assembled through a hierarchical process, with 
small structures (M ~ 10^ Mq; about 10^^-10^^ of the present galactic 
mass) forming first. Such structures are estimated to be a few hundred 
pc in size, and to have virial velocities of order 10 km/s. If magnetic 
fields grow on the scales of the largest eddies in these systems, the 
result would not be very different from eddy scale fields in a massive 
galactic halo. But the problem of eliminating the small scale fields 
and growing large scale fields remains. 

In summary, amplification of the magnetic field in the protogalaxy 
could address problems with the disk dynamo related to its operation 
at very low field strength. The flux expulsion problem, however, still 
remains even at the present field strength. 

13 Generation of small scale magnetic 
fields by turbulence 

Earlier, in §6, we discussed how, based on turbulence with helicity 
flows, the a — Q dynamo can generate magnetic fields with net flux, in- 
side a disc. In the last section we considered how magnetic fields could 
be produced by dynamo action in turbulence with no helical flows. An 
astrophysical application of such a dynamo raises the possibility of 
producing substantial pregalactic fields in the protogalaxy. Another 
application is the production of fields in the intracluster plasma in 
clusters of galaxies. (We do not consider the latter, as the main point 
of this review is to uncover the processes involved in generating cosmic 
magnetic fields, and in this review we have reduced the problem to 
the challenge of producing a magnetic field in our Galaxy.) 

In the case of the protogalaxy, we have seen that magnetic fields 
produced by nonhelical turbulence have no mean flux overall, but 
could have mean flux in large subregions. Other subregions may have 
a mean flux with a different direction. However, such a field can be 
the source of the Galactic field which itself has no overall net flux; 
only mean flux over large regions. Magnetic fields in protogalaxies are 
generated by random, nonhelical, small-scale turbulence. 

Such a theory of field generation has been developed by a number of 
people (Kraichnan and Nagarajan 1967, Kazantsev, 1968, Kulsrud Sz 
Anderson 1992, Subramanian & Barrow, 2002, Boldyrev &; Cattaneo 
2004). Here we review the Kulsrud-Anderson theory in some detail. 

Assume that the statistics of the turbulence is given by a turbulent 
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spectrum such as that of Kolmogorov. For simpHcity, assume that the 
turbulence is homogeneous, isotropic, and steady in time. Assume 
that initially there is some very weak magnetic field, whose statistics 
are not necessarily isotropic, but are homogeneous. 

Then one can evolve the magnetic field by integrating the magnetic 
differential equation (2), with the resistive term set to zero, (In almost 
all cosmic situations it is permissible to neglect the resistivity except at 
extremely small scales) . The velocity in this equation can be assumed 
to be given at least statistically and unaffected by the magnetic forces. 
This is the kinematic assumption. 

Consider the system to be enclosed in a large periodic box so that 
both the velocity and the magnetic field can be described by Fourier 
transforms 

V = y d^kvi,it)e'^-^ , 

B = y d^kBk(i)e*^ (35) 

Since the turbulent velocities are clearly random, we describe them 
statistically by an ensemble average, and make use of the random 
phase approximation (Sagdeev and Galeev, 1969, Kulsrud 2005) 

< Vk/(t)vk(t') >= J{k){l - kk)(5(k' - k)Sr{t - t'). (36) 

Here the angular brackets denote the ensemble average, of the veloc- 
ities at different wave numbers and times, I is the unit dyadic, and 
the hat denotes a unit vector. The delta function in k' — k results 
from the random phase approximation. Its significance is that there 
are a large number of independently positioned velocity waves in the 
large box and taking the Fourier transform of these waves yields a 
delta function factor between different wave numbers. S-j-it ~ t') is the 
correlation function of the velocities at a fixed position at different 
times. If we normalize this function to unity when t = t', then J{k) 
is the energy spectrum of the turbulence 

1 < v2 >= J d^kJ{k). (37) 

As a result of the turbulent velocities and the magnetic differential 
equation the one dimensional magnetic energy spectrum M{k), given 

by 

— = / M{k)dk (38) 

Btt J 
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evolves in time according to the one dimensional mode-coupling equa- 
tion 

dM{k) 
dt 

The kernel K{k, ko) is 



= J K{k,ko)M{ko)dkQ-2pk'^M{k). (39) 



K{k,ko) = 2TTTk^ f ^-^^^-^(k^ + kl - kkocose)j{ki), (40) 
J ki 

where ki is defined by 

kl = k^ + kl-2kkQCOse. (41) 

The quantity (3 is the same turbulent diffusion coefficient as in the 
a — theory and is 

P = It I J(fc")dV. (42) 



3 

r is the decorrelation time given by 



T 



POO 

/ 6r{t)dt. (43) 
Jo 



(The mode coupling equation is strictly valid only if r is sufficiently 
small compared to the rate of evolution of the magnetic spectrum. 
This is called the short correlation-time approximation.) 

The derivation of the appropriate form of the mode coupling equa- 
tion is given in a number of places (Kraichnan & Nagarajan 1967, 
Kazantsev, 1968, Kulsrud &; Anderson 1992). The form given here 
was derived in the paper of Kulsrud &; Anderson. In this derivation, 
no assumption as to the isotropy of the magnetic field statistics is 
made. M(k) refers to the three dimensional magnetic spectrum aver- 
aged over all solid angles. 

This equation has a number of important implications. The most 
important one relates to the growth of the total magnetic energy 

S = J M{k)dk. (44) 

If one integrates equation (39) over k one finds the evolution equa- 
tion for £, 

f = 2.£, (45) 
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where 

1=^J k'^J{k)d^\i =^ J k^J{k)dk. (46) 

A very similar equation was discovered by Batchelor (1950). 

This equation has an easy interpretation for Kolmogorov turbu- 
lence. Namely, eddy velocities with wave number k have a turn over 
rate of kv where = J(fc)d^k is the mean square turbulent velocity in 
the range d^k. The decorrelation time r is ~ l/fev. The change 5B/B 
during the decorrelation time, is plus or minus the rotation angle of 
the eddy M = kvr. Thus, < {dB/Bf/r >= k'^v'^r = k'^J{k)Td^k, 
the contribution of the integral in d^k to 7 in equation (46). In short, 
every turnover of any eddy produces a relative change in the mag- 
netic energy by a factor of order two. Note also, that each eddy acts 
independently of the others in increasing the magnetic energy. 

How does this change in magnetic energy produced by a single 
eddy manifest itself in the change of the actual field? Surprisingly, 
the magnetic field change produced by any eddy occurs on scales much 
smaller than the eddy itself. It had been previously assumed that the 
magnetic energy is produced on the same scale as the eddy. (See for 
example the treatment in RSS, 1988). This is because, in the presence 
of a uniform or large scale eddy, the main interaction was previously 
assumed to be between the eddy and this large scale field. Actu- 
ally the stronger interaction is between the eddy and the small scale 
magnetic field, and this produces energy at twice the magnetic wave 
number. As a result, magnetic energy rapidly propagates to larger and 
larger wave numbers. The wave number, at which the amplified en- 
ergy peaks, increases exponentially, and the field develops a very much 
finer structure than that of the eddy. This result might be anticipated 
since by flux freezing the magnetic field follows the displacement of 
the plasma rather than its velocity. Since the displacement is the inte- 
gral of the velocity one expects the displacement to also exponentially 
develop fine structure. 

These results are demonstrated analytically in Kulsrud & Ander- 
son (1992) as follows: Consider that the turbulent — eddy scales are 
larger than 1/fci. Then expand the mode coupling equation (39) for 
larger values of the magnetic wave numbers, k and /cq, where we see 
from the notation in the mode coupling equation (39) that magnetic 
energy at the large wave number ko is transferred to the also large 
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wave number k. The result is the simple partial differential equation 
'^-lf*^5¥-2*l^ + eMV2*%M (47) 



dt 5 V dk"^ dk 

where A5 = rjsc^/^T^ is the magnetic diffusivity and rjs is the Spitzer 
resistivity. (The resistivity term is added by hand.) 

Without the resistivity term, this equation has a Green's function 
solution 

M{k,t)=f M{-^,t')G(^,t-t'))dt', (48) 
where G is the Green's function 

This solution expresses M{k, t) for general k in terms of its time vari- 
ation at some reference wave number kj-ef- 

Thus, if M{k/kref,t') is a delta function in time then as t increases 
the magnetic energy increases as k^^"^ up to a peak near e^/^, beyond 
which it decreases. At the same time it increases at a fixed k by the 
exponential factor e^^^^^. Thus, the energy grows exponentially at 
fixed k and spreads out over a range that also increases exponentially. 
The combination of these two effects, local growth and growing range, 
leads to the growth in the total energy at the rate 27. 

Since the scale to which the magnetic energy extends gets small 
so fast, one expects it to reach a resistive scale, i.e. a scale at which 
the resistivity becomes important. Restoring the resistivity term in 
equation (47), one finds a purely growing solution 



M{k,t) = const e(3/4)7*A;3/2^Q (v 
where Kq is the Bessel function of the second kind. This shows that the 



wave number corresponding to the resistive scale is krj ~ ^/^y/WXs■ At 
larger wavelengths M drops off exponentially as the magnetic energy 
is destroyed by resistivity. However, even with resistive damping the 
total magnetic energy still increases exponentially, but at the slower 
rate 87/4. 

One can fix the constant so that the energy at some k,.(,f is 
M{kref,t) = e''^^!^ . If one does that, then one finds that at the time 
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the small scale energy first reaches A;^, the total magnetic energy is 
larger than that at kref, krefM{kref), by the factor v/kref^s-, which 
is the effective magnetic Reynold's number. 

For application to protogalactic magnetic fields, it is of interest to 
estimate the order of magnitude of the growth rate of the small scale 
turbulent energy while still in the kinematic limit. 

In this case, the growth rate is the rate of turnover of the smallest 
eddy, which we take to have wave number ki,. The eddy turnover time 
is {kvk)~^ where Vk is the velocity at the eddy with wave number k. For 
Kolmogorov turbulence ~ k~^/^ so the turnover rate is proportional 
to A;^/^. The ratio of the turnover rate of the smallest eddy to that 
of the largest eddy is (k/kof^^ where k^ is the wave number of the 
largest eddy. Now, if the turbulent energy is in equipartition with 
the gravitational energy, it is easy to see that the turnover rate of the 
largest eddy is equal to the free fall time of the protogalaxy. (This 
is also equal to the decay time of the turbulence.) But /c^/fco = R^^^ 
where R = 2iTVQ/kQV is the fluid Reynold's number, v k. Vfhi is the 
kinematic viscosity for an ion thermal speed Vf^ and mean free path 
£ = 1/na, where n is the ion density and a is the Coulomb cross 
section. 

Combining these results and setting the free fall time equal to td 
we get 

kov 

or 

/ 27r?;o \ vqLo L r— 

^^-=(l^J=^"l = ^^' 

where Lq = lir jk^ and we take vq w Vth from equipartition. 

Now more specifically we take the baryonic mass of the protogalaxy 
to be Mb = 10^^ solar masses and the protogalactic radius to be one 
hundred kiloparsecs. The dark mass may be ten times this. Then the 
Jeans temperature Tj Ri GMmn /Lis 10^ degrees Kelvin or Tg^ ~ 10^ 
eV. The coulomb cross section a ^ 10~^^/r^ 10~^^cm^. 

The baryonic density is given by 47rL^/3n = Mb orn ~ 10 cm 
Therefore from our equation we get 

2'yTD « 55. (52) 

( A more precise estimate and for a general galaxy yields 2^td ~ 
70/ y/Mii and depends only on the baryonic mass. Mn is the baryonic 
mass of the protogalaxy in units of 10^^ solar masses.) 
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Thus, we see that the smahest eddy turns over seventy times during 
the free fall time and in the kinematic limit the small scale magnetic 
energy grows by a factor 2 x lO^'^ while the rms field strength would 
increase by 1.7 x 10^^. By equation (33), the Biermann battery gener- 
ates a field of roughly 10~^^ gauss on this scale, which leads to a field 
of 3 X lO""^ gauss. However, the energy density of the smallest eddy is 
only \/2{nmHVQ/ R}^"^) ~ 10~^^ergs /cw? so the kinematic assump- 
tion breaks down and saturation of both the smallest eddy and the 
magnetic field must occur. 

This amplification proceeds even faster in the smaller, cooler, less 
viscous structure such as the ~ IQ^Mq halos, with which galaxy for- 
mation is generally believed to begin. However, the estimate of the 7 
is very dependent on the turbulence theory which is highly uncertain. 
The estimate given in Kulsrud et al (1997) gives a value of which 
is larger by 300 and this would produce an increase of the field energy 
by 6 X 10^^, and a resulting field strength of about one gauss, and a 
much larger magnetic energy density than the smallest eddy. Again, 
saturation of the magnetic field must occur. Saturation is discussed 
in the next section. 

14 The saturation of the small scale 
magnetic fields 

In the previous section we discussed magnetic fields generated by tur- 
bulence and found that the small scale magnetic energy, that at scales 
below the smallest turbulent eddy, increases very rapidly, doubling 
essentially every turnover time of this smallest eddy. In the kinematic 
limit, the energy propagates down to the resistive scale after which 
time it is larger than the magnetic energy at the eddy scale by a 
factor of the magnetic Reynolds number. 

But such a large energy is clearly beyond the kinematic limit. The 
magnetic field must produce nonnegligible velocities which saturate 
the magnetic energy at a much smaller value. For large magnetic 
Reynolds numbers the kinematic limit breaks down before the resistive 
scale is reached. This is because the small scale magnetic energy 
saturates. 

This effect was pointed out in a different way by Cattaneo and 
Vainshtein (1992), and Diamond & Gruzinov (1997). (Also see Gruzi- 
nov h Diamond (1994, 1995). The statement that these authors make 
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is that the a coefficients of the a — Vt theory are suppressed when the 
large scale magnetic energy is larger than the turbulent kinetic energy 
divided by magnetic Reynolds number. This puts a severe restric- 
tion on any dynamo theory since the magnetic Reynolds number for 
galactic turbulence can be as large as 10^^. However, if we view their 
statement from the point of view of small scale theory and replace 
M{kref) by the large scale field taken at the eddy scale, we see that if 
their condition is satisfied and M{kref) is smaller than the eddy en- 
ergy divided by the magnetic Reynold's number, then the small scale 
fields reach the resistivity scale before saturation. Otherwise, they are 
not able to reach the resistive scale. Instead, energy equipartition of 
the total magnetic energy of the small scale fields is reached with the 
smallest eddy, and its dynamo effect is modified. But after equiparti- 
tion is reached, the small scale energy stops growing, and never gets 
much larger than the eddy energy. Thus, the suppression of the dy- 
namo coefficients is not limited by the large factor of the magnetic 
Reynold's number. 

To understand the saturation process it is necessary to understand 
the structure of the small scale fields. This structure has been found 
numerically by Cowley et al. (2004), motivated by a theoretical hy- 
pothesis of Cowley. They found that the magnetic fields on small 
scales are far from isotropic, and indeed consist of folded magnetic 
field loops whose length is that of the eddy size while the thickness of 
the folds is very small. (See Figure 3 where such folds are displayed.) 
The folded nature does not emerge from the mode coupling equation 
(39) of the last section, since this equation only refers to M{k), the 
magnetic energy at a given k averaged over all solid angles. This av- 
eraging hides the structure of the fields. As can be inferred from the 
figure, the Fourier transform of a folded field assigns a wave number 
equal to the reciprocal thickness of the fold and the long length drops 
out. 



On the basis of these folded fields it easy to understand how mag- 
netic fields on scales much smaller than the inner turbulent scale can 
be amplified. To the folded fields on such small scales the smallest 
turbulent eddy appears quite large and the magnetic folds only feel 
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Figure 3: A folded field. It is amplified by the stretching velocities v's, due 
to the smallest turbulent eddy. At equipartition it can be unfolded by the 
velocities u's at the Alfven speed if Braginski viscosity holds. 

the effect of a part of the eddy, the shearing part. The shear of this 
part appears constant in space. This shear can also be envisioned 
as an incompressible flow that is compressing in one direction and 
expanding in a perpendicular direction. If the compression direction 
is perpendicular to the fold, the fold will be squeezed and stretched, 
and its magnetic field strength will be increased. Of course, if the 
compression is parallel to the fold and the stretching is perpendicular, 
then the field strength will be decreased. But for random shifts in 
direction of the shearing flow, the net effect will be an increase in the 
mean square of the field strength of the folds as well as a decrease in 
their thickness. This mechanism exactly mirrors the more complex 
mode-coupling calculations of the last section. The validity of this 
modeling of the amplification is borne out by numerical simulations. 
(Maron et al 2004, Schekochihin & Cowley 2004). 

Prom inspection of Figure 3 one can also see how the folded mag- 
netic fields should tend to saturate. Indeed, neglecting any viscosity, 
we see that the folded fields have tension at their ends and this tension 
should produce a shortening of the fields at a velocity u equal to the 
effective Alfven speed based on the strength B of the folded fields, 
divided by their length. To see this, note that the current density j 
perpendicular to the fold is of order B / 4:Tt5 where 5 is the thickness of 
a single fold. The field perpendicular to the fold is also not zero as ap- 
pears in the figure, but is of order {S/L)B where L is the length of the 
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fold. Thus, the force on the plasma along the fold is F = {B)'^/4TrL, 
and this leads to a velocity u ~ {F/ p) x t where t ~ L/u is the time to 
unfold. Thus, v? « {B)^ /A-Kp and so u is approximately equal to the 
Alfven speed. This implies that at saturation, when the small scale 
magnetic energy is the same as the kinetic energy of the eddy, the 
unfolding time is comparable to the eddy turnover time, and to the 
exponentiation time of the field; the folds no longer grow or decay. In 
other words, the small scale fields saturate by unfolding as fast as they 
grow when their energy is in equipartition with the kinetic energy of 
the smallest eddy. 

The same level of saturation is found in the numerical simulations 
but the saturation seems to have a different origin since the folded 
fields continue to elongate in these simulations. The different behav- 
ior is attributable to the fact that the numerical simulations include 
fluid viscosity and this prevents the unwinding motions since the un- 
winding velocity on one fold is in the opposite direction to that on the 
next fold. Because of the very great thinness of the folds, the viscous 
forces between these flows is enormous. (Recall that viscosity termi- 
nates the turbulent spectrum at the scale of the smallest kinetic eddy 
whose size is much larger than the thickness of the folds.) However, 
the viscosity used in the simulations is ordinary scalar hydrodynamic 
viscosity. Because of the magnetic field, the folds can not transmit 
shear forces across the magnetic field for distances larger than the ion 
gyration radius, a distance generally much thinner than the thickness 
of the folds. The proper viscosity to use in this case is the Braginski 
tensor viscosity, which docs not allow such a perpendicular transfer of 
shear forces. (Braginski 1965, Malyshkin & Kulsrud 2002.) 

These pictures for the evolution of the folded fields, the one based 
on ordinary viscosity and the other based on Braginski viscosity, lead 
to the same equipartition of the fields generated by the smallest eddy. 
However, they should lead to differences when the next larger eddy 
starts to play a role. For ordinary viscosity the next larger eddy 
should stretch the length of the folds to equal itself. For Braginski 
viscosity this tendency should be suppressed and the magnetic field 
should develop a more isotropic structure on the scale of the next 
eddy. 

In this case, the resulting field should have a much larger scale, 
possibly large enough to appear very coherent. In the case of the pro- 
togalaxy, the resulting field supplied to the galactic disc should have a 
scale comparable with observations of the present galactic fields. On 
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the other hand, for normal viscosity the magnetic field should end up 
very thin and with folds of the length of the largest eddy. Its energy 
should be in equipartition with the turbulent energy up to the eddy 
for which nonlinear growth can arise. This energy can be a fair frac- 
tion of the turbulent energy and thus of the binding energy of the 
entire protogalaxy, so it is possible that the magnetic field will alter 
the collapse to the disc. 

At the moment numerical simulations of magnetic field evolution 
in the presence of turbulence employ ordinary viscosity because Bra- 
ginski viscosity is numerically difficult to handle. As yet, the results 
from Braginski viscosity are all based on theoretical estimations, and 
we must wait for more sophisticated numerical programs, yet to be de- 
veloped, before the saturation is properly understood. However, it is 
generally agreed that Braginski viscosity is the correct viscosity when 
the ion gyration radius is very small compared to all the other scales. 

An additional problem has arisen because the Braginski viscosity 
in very weak fields leads to instabilities. (Schekochihin & Cowley 2004, 
2006; Schekochihin et. al. 2005a, b.) The instabilities arise because in 
very weak fields the changing magnetic fields produce anisotropies in 
the pressure tensor that drive fire hose and mirror instabilities. These 
are normally suppressed in stronger fields. How these instabilities 
saturate is still under discussion. 

Assuming the folded fields saturate by unfolding, Kulsrud has at- 
tempted to build an intuitive model for the evolution of the proto- 
galactic magnetic fields in Kolmogorov turbulence. (Kulsrud 2005) 
Assuming reasonable parameters and top down formation of the pro- 
togalaxy, he finds that the fields are progressively generated by larger 
and larger eddies, and in the free fall time the size of the most inten- 
sive fields can reach about one thirtieth the size of the entire proto- 
galaxy. Also, the field strength, after compression into the disc, is in 
the microgauss range. This speculation does not take into account the 
instability problem. 

Thus, the saturated fate of magnetic fields generated by turbu- 
lence, must wait till numerical simulations with the appropriate vis- 
cosity can be developed and the problems with the nonlinear Braginski 
instabilities understood. 

However, because the validity of the a — J7 theory is equally in 
doubt, it is of interest to consider the implications of a protogalactic 
field, and its further evolution after the galactic disc has formed. 
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15 History of the evolution of a pri- 
mordial magnetic field 

In the last two sections we consider the possibility that the proto- 
galaxy could generate a strong enough magnetic field that it might be 
considered a pregalactic field. What would be the evolution of a pre- 
galactic field which existed prior to the formation of the Galactic disc, 
and is subject to the astrophysics of the Galactic disc? How would it 
survive the criticisms of Parker discussed in §10? 

The answer to this question was addressed in a paper by Howard 
& Kulsrud, (1997). A condensed summary of this theory is given by 
Kulsrud in Cosmic Magnetic Fields edited by Wielebinski and Beck 
(Kulsrud, 2006). Here an even briefer treatment of this paper is given. 



First, imagine a uniform magnetic field present in the sphere of 
plasma that is to collapse to form the Galactic disc. The behavior of 
the field is sketched in Figure 4. In panel (a) the uniform field is dis- 
played in the protogalactic sphere after it has separated from Hubble 
expansion. In panel (b) the plasma collapses to a smaller virialized 
radius dragging the field lines with it. However, it is seen that they 
are still connected to the surrounding intergalactic medium. Possi- 
ble radii for these spheres might be 100 kpc and 10 kpc. If the field 
makes a finite angle to the rotation axis, the horizontal component of 
the magnetic field would be amplified as the inverse radius squared, 
or 10^. Next, in panel (c), the second sphere has collapsed to a disc 
which might be 100 pc thick. In this case the horizontal component 
would be amplified again by another factor of 100, and the final field 
strength would be increased by 10^ over that of the original primordial 
field. The vertical component would be changed by only a factor of 
10^. 

The resulting field lines are displayed in this third panel. It is seen 
that some lines, of type A, enter the disc from above and also leave in 
the above direction. However, some others, type B, enter from below 
and leave above. Lines of type C both enter and leave below. Lines of 
type B are embedded in the disc and cannot escape without violating 
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Figure 4: The early evolution of a primordial field, (a) The uniform field 
in the protogalaxy before collapse, (b) The field lines after the protogalaxy 
has shrunk spherically and they have been drawn inward from the extra- 
galactic region. (c)The field lines after the disc forms, the lines A entering 
and leaving from above, lines B entering from below and leaving above and 
lines C entering and leaving from below, (d) The solid line before ambipolar 
diffusion and the dotted hne after. 
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flux freezing. The other hnes do leave by ambipolar diffusion. 

The subsequent behavior depends on the structure of the inter- 
stehar medium. It is generally supposed that the magnetic field is 
embedded in many clouds. It threads its way through these clouds 
and eventually finds its way out of the disc as in Figure 5. The clouds 
hold the magnetic field in the disk balancing its outward force by the 
gravitational field of the stars acting on the mass of the clouds. (In 
addition, there is also a cosmic ray pressure which is held in by the 
clouds through the action of the magnetic field and indirectly adds to 
the outward pressure exerted by the magnetic field.) The low density 
material between these clouds is assumed to be very light and unable 
to hold the field and cosmic ray pressure down, so the field lines bow 
up between the clouds. 

Thus, the force balance is between the outward pressure of the 
magnetic field and cosmic rays acting everywhere including the space 
between the clouds. This is balanced by the gravitational force acting 
only on the clouds so that the effect of the field and cosmic rays is 
concentrated and amplified in its action on the clouds. If the filling 
factor of the clouds in the interstellar medium is /, then the force on 
the clouds is intensified by a factor of 1//. The effect of the magnetic 
field is further increased by the fact that the clouds are only very par- 
tially ionized, generally by 10^^ or 10~^. There is no direct magnetic 
force on the neutrals, but it is transferred to them by the collisions 
between the ions and neutrals. This force only acts on the ions, and 
the field lines locked on them, slide through the neutrals, and through 
the clouds, at some ambipolar velocity. Because of the amplification 
factors it is possible for the ions and the field lines embedded in them 
to slip entirely across a cloud. 

The situation is not entirely simple since the time to slip through a 
cloud is much longer than the lifetime of the cloud. (The slippage time 
is of the order of a billion years while the lifetime of a cloud before 
it hits another cloud and is destroyed, is in the range of ten million 
years.) However, this complication is ignored since after a cloud is 
destroyed its plasma is turned into a rarefied plasma in which the 
field is still embedded by flux freezing. After a short time a new cloud 
reforms out of this plasma drawing the field lines back into itself with 
the ionization dropping back to its very small value again. After the 
lifetime of a single cloud the matter has worked its way outward from 
the midplane of the disc, in a non random way. For a given line the 
amount of plasma that is below it (relative to the disc) has increased. 



67 



B 




Figure 5: A line of force penetrates many clouds and finally reaches a last 
cloud through it diffuses and leaves the disc. 



The amount of increase is given by the amount the ionized matter, 
holding down the field line moves in order to balance the magnetic 
and cosmic ray forces against the ion neutral collision rate. 



Finally, the field line manages to slip through the last cloud it 
threads before it leaves the disc, as is shown in Figure 5. Now, this 
motion will remove lines of type A and C, shown in panel (c) of figure 
4, out of the disc. But it cannot move field lines of type B out the disc 
because of their topology. However, when lines of type B slip through 
the last cloud mentioned above, the result is that they shorten their 
horizontal extent in the disc as seen in panel (d) of Figure 4. 

In discussing the behavior of magnetic fields, the authors ignore 
any interstellar turbulence or dynamo action. This enables them us 
to isolate and simplify the ambipolar effects on the magnetic field and 
its lines of force. 

The only motions of the interstellar medium that they include are 
the galactic differential rotation and the vertical ambipolar velocity. 
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In cylindrical coordinates these velocities are 



V = m{R)e + Vzi, 



Vz = -K- 



dZ 



K^^lliM^ (53) 

Vz is the ambipolar velocity needed to transfer the cosmic ray and 
magnetic forces to the neutrals in the cloud, pi is the ion density, 
is the ion- neutral collision frequency and the factor P/a is the ratio 
of the cosmic ray pressure to the magnetic pressure, which is assumed 
constant with height. 

To see how the fields evolve under these motions ignore the vertical 
component of the magnetic field Bz- Then, with these velocities, the 
magnetic differential equations for the field components in a frame 
following the galactic rotation are (Kulsrud 1990), 

dBR_dBn dBR_ d 
dBe dBe ^ dBe ^ ,t/ r ^ ^ r ika^ 

Let the half thickness of the disc be D. Outside of the disc it 
is now appropriate to take these components to vanish since they are 
embedded in the very tenuous plasma that is ionized in the clouds and 
has negligible mass. (The weakening of the horizontal components, 
Br and Bg of the field by the Vz motion is due to the shortening of 
the field line length L, where the horizontal field strength is given by 
B = (L/D)Bz, and Bz is constant.) 

To form an idea as to the solution assume that the fields are 
parabolic in Z. The horizontal dependence is not given since a single 
column of plasma is followed as it is carried around the disc by the 
differential rotation. 

Thus, take 

BR{Z,t)=BR{0,t){l~Z^/D^), 
Be{Z, t) = Be{0, t){l - Z^D^), (55) 

and solve for the time dependence of the fields at Z = 0, the coeffi- 
cients of the parabolic terms on the right hand side of this equation 
(55). Taking advantage of the fact that shear quickly amplifies Bg 



69 



relative to Br, the solution at the midplane Z = can be written 
with sufficient accuracy as 

Bn 



[l + (2F^i/3D)J72i3]i/2' 

-Bint 

Be = n^, (56) 

[1 + {2VDi/3D)nH^f^ 

where Bi is the initial value of Bji, Vdi = KBf/2TTD is the initial 
value of Vd, and assume that initially Bg = 0. Bg is plotted versus 
time in Figure 6 for several values of Bi. 



These equations display the behavior of the field in the midplane 
of the disc following a rotating fluid element. Their behavior can be 
described qualitatively as follows. 

Consider some fluid element that initially has a relatively weak 
initial B^ component Bi. At first the only action is the differential 
stretching that leads to a linear growth of the toroidal field Bg. Then, 
as Bg grows stronger, the ambipolar velocity starts to grow and be- 
come important. This leads to a weakening of Br, but because of 
differential stretching, Bg continues to grow for a time. Eventually, 
Bji becomes too weak for stretching to overcome the weakening of the 
Bg, and it too starts to decrease. Eventually, when the second terms 
in the brackets in the denominator dominate, the two components 
decrease as powers of t. 

The stronger toroidal component then varies as 



It is easy to see that this field is just such that the ambipolar dif- 
fusion velocity, which is proportional to B^, moves the field a distance 
D in the time t. From the definition of v^i one sees that at any later 
time, t, which is sufficiently long, the field strength is independent of 
its initial value and only depends on the initial sign of Br. Thus, the 
entire magnetic field throughout the disc becomes a constant in space 
except for its sign. 
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Figure 6: The evolution of the field strength under differential rotation and 
ambipolar diffusion for initial field strengths of Bl =,0.1,0.3 and 1 micro- 
gauss. The ordinate is B in microgauss and the abscissa is time in biUions 
of years. 
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One can take this sign into account by noting that, according to 
the model for the primordial field, Bji = Bi cos 9 and at a later time 
increases by Q,{R)t. Because of the variation of Q, with R the magnetic 
field at the time t and fixed 9 has a square wave behavior with R. The 
spatial period is Ai? = {(Kl/dR)~^. [See equation (25).] This should 
be roughly constant with R near the sun and is approximately 100 
parsecs. This is the wound up field discussed in §10, in the context of 
Parker's discussion of the winding up. 

Such a field cannot agree with the pulsar rotation measures because 
most of the Faraday rotation in the intervening interstellar medium 
would cancel out. However, assume that the initial magnetic field 
is not uniform, but had a gradient with a scale comparable to the 
protogalaxy. Then after collapse Br is not sinusoidal. Bg gets its sign 
from the initial sign of Br in the same piece of plasma. Thus, the 
resultant radial variation of the final field Bq is square wave with the 
radial extent of one sign being different from the radial extent of the 
other sign. 

Now, for such a field the Faraday rotations do not average out. 
Further, a field with such variations would change on too small a scale 
to be resolved by the analyses of rotation measures available at the 
present time. This is just the field introduced in §10 to counter the 
objection to the winding up of the pregalactic field. 

It must be remembered that in this model all dynamo actions and 
dynamics are ignored. Including these would lead to a different field 
structure. 

What should be taken seriously is that the suggested pregalac- 
tic field could give the dynamo a proper start to overcome the fiux 
loss problem. The initial pregalactic field could be relatively weak, 

and then be developed to sufficient strength by the standard a O 

dynamo to satisfy the necessary vacuum conditions as discussed in §9. 

The purpose of the work in Howard &: Kulsrud is not to directly 
prove that the magnetic field is primordial but to examine the con- 
sequences of its being so. It potentially could give a strong field but 
would not give the Grand Design shape that generally emerges from 

and a Q theory. However, it cannot be so easily ruled out on these 

grounds because as shown in §10 the perturbations produced by com- 
pression of the field in the spiral arms can also give a similar pattern 
to the field lines. 
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16 Ext r agalactia magnetic fields 



Most of this review is concerned with the magnetic field of our Galaxy 
and its possible origin. The problems encountered in explaining the 
origin of the Galactic field are not universal - galaxies with higher 
star formation rates or smaller gravitational fields could expel flux - 
but these are not representative of all galaxies. 

A much more difficult problem is the origin of the magnetic fields 
observed in the intracluster medium (ICM) of clusters of galaxies (Car- 
illi and Taylor 2002, Ensslin et al. 2005). These fields are coherent on 
scales of several kpc and in the cores of clusters appear to be several 
microgauss in strength, comparable to that of galactic fields. It can be 
inferred from the enrichment of the ICM with heavy elements, and the 
paucity of gas rich galaxies in clusters, that interstellar gas is stripped 
from galaxies and mixed into the ICM. However, the gas density of 
the ICM is 2 - 3 orders of magnitude lower than the mean density in 
galaxies, so the field in the ICM must be amplified by some mecha- 
nism to bring it back up to microgauss strength. The large intrinsic 
scales and lack of defined rotation of clusters make the operation of 
a mean field dynamo problematic. Various other means to produce 
the cluster field are discussed in Carilli & Taylor (2002), Schekochihin 
et al (2005a,b), and Ensslin et al. (2005), but none intrinsic to the 
cluster seem promising at this time. 

However, a proposal that seems promising is based on the observa- 
tion that a large magnetic flux fills each of the giant radio jets whose 
volume is enormous. If when such jets disband, their fields could be 
dispersed uniformly throughout intergalactic space, then there would 
be enough flux to supply significant pregalactic fields to every galaxy. 
(Kronberg, Duften & Colgate 2001). 

This was first pointed out by Daly & Loeb (1990). This result is 
based on the assumption that radio galaxies are rather short lived, 
and new ones emerge every few hundred million years or so, each new 
one producing more fiux. 

When the radio jets terminate their life their plasma expands into 
the intergalactic medium taking their flux with it. If the intergalac- 
tic medium is turbulent, this flux then merges with the surround- 
ing intergalactic plasma mixing its flux with it. This intergalactic 
plasma subsequently collapses to produce galaxies which are thus al- 
ready magnetized when they form. That the total flux in radio jets 
is adequate to produce the observed flux in clusters and galaxies, is a 



73 



striking numerical coincidence. How this happens, the expansion and 
mixing of the flux, is not yet understood. However, given the problems 
with other theories for the origin of cosmic fields, this origin must be 
taken seriously and developed further. 

Much work has already been done on this problem by Colgate &; 
Li (2000) and others. They have developed a preliminary theory to 
explain the origin of the magnetic fields of the jets. Their theory 
is based on ideas developed by plasma physicists to explain some of 
the interesting fusion experiments. The principle mechanism which 
is employed is Taylor relaxation (Taylor 1974, 1986), in which, by 
rapid reconnection, a spherical plasma relaxes to a particularly sta- 
ble configuration by converting toroidal flux to poloidal flux and vice 
versa. Determining whether this is possible on the large megaparsec 
scale of radio jets seems to depend on increasing our understanding 
of the fast magnetic reconnection process. Study of this is a field now 
under intense development. For reviews see Priest &; Forbes (2000), 
Biskamp (2000), Kronberg et al. 2001, Furlanetto & Loeb 2001, and 
Gopal-Krishna & Witta 2003). 

Presumably, the fields in the jets originate in accretion disks. Many 
of the processes already discussed in this article - the Biermann bat- 
tery and , various types of dynamo - should be able to generate and 
amplify magnetic fields in disks much more rapidly than in galaxies, 
particularly if the disks are hot and collisionless (Krolik &; Zweibel 
2006). Once these fields are ejected into the jets, the problem remains 
to amplify them, and generate some semblance of a coherent field. 

17 Summary and conclusions 

The origin of cosmic magnetic fields breaks up into three stages: (1) 
The generation of very weak seed fields from an absolutely zero mag- 
netic field. (2) The amplification of these seed fields to dynamically 
interesting fields. (3) The final alteration and further amplification to 
the fields we observe today in galaxies and clusters of galaxies. 

The first stage is generally thought to occur in the primordial era, 
before and during the time when galaxies and their discs form. There 
are a variety of mechanisms for their production, the Biermann bat- 
tery mechanism being the most popular. Another mechanism is the 
generation of fields in stars (by the Biermann battery) and the later 
expulsion of their flux into the interstellar medium. In short, there 
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seems to be no difficulty witli finding mechanisms to produce the first 
stage seed fields, and not much debate about which is the most im- 
portant. It is just generally accepted that they are present as needed. 

The main debate centers on the second stage, as to when do dy- 
namos start to amplify seed fields. Does this occur pregalactically, 
before galactic discs are formed, or does it occur over the life time 
of galaxies? The bulk of opinion during the past decades favors the 

latter choice. It has been generally accepted that the a f2 dynamo 

takes these seed fields and amplifies them up to the currently observed 
fields. There are several difficulties with this theory. These difficul- 
ties are primarily associated with the flux freezing condition that says 
that the total flux cannot change from zero to a finite value in fixed 
amount of plasma. This condition is satisfied in standard disc-dynamo 
theory by imposing boundary conditions that imply that during any 
amplification of the field by a finite amount, a finite amount of flux is 
removed from the disc to infinity, leaving behind an increased amount 
of fiux of the opposite sign. Thus flux is conserved in toto but not in 
the disc region alone. 

The boundary conditions in the a theory that enforce this 

condition are the vacuum boundary conditions. These conditions are 
formally derived by assuming that turbulent diffusion outside of the 
disc is much larger than that in the disc. This outside turbulent 
diffusion is based on observations of velocities and correlation lengths 
made in the halo. But these measurements are made when the halo 
has very little matter in it. If these measurements were made at a time 
when tubes of force carrying matter with them were diffusing out of 
the disc, the gravitational force on the matter would strongly modify 
these velocities. In fact, for very weak fields, the physics of the removal 
is difficult to justify, by models which treat the interstellar medium 
as horizontally homogeneous. The conclusion that we can draw is 

that the a theory that treats the the interstellar medium as 

horizontally homogeneous is astrophysically incorrect and this is likely 
to be true even in the case of strong fields. 

It is conceivable that a more sophisticated treatment that takes 
into account the non homogeneity of the interstellar medium might 
be justifiable. In this case, when the magnetic field fiux leaves, the 
field lines and matter would form arcs in which the matter slides down 
the tubes of force releasing the remainder of the tube to more freely 
diffuse away. But this is clearly not possible when the field is very 
weak, so weak that its forces do not affect the motion of the matter. 
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which would be unaware of the direction to shde down. 

Thus, the standard theory is probably not able to amplify these 
weak fields, and some other mechanism is necessary to provide this 
intermediate step in the origin of cosmic fields. That is to say, a pre- 
galactic mechanism might be necessary to provide fields of reasonable 

strength for the non homogeneous a disc dynamo to act on. In 

this sense with a pregalactic origin to get started such an a 0, 

dynamo could play a crucial role in the magnetic field origin. 

Even if the flux removal problem is solved, the a Q dynamo is 

not necessarily viable. There remains a gap between the calculations, 
both analytical and numerical, which show tremendous growth of fields 

at the resistive scale, and the mean field a 17 dynamo models, which 

involve the transport of magnetic flux on the much larger supernova 
remnant scales. It is supposed that the small scale fields can saturate, 
leading to dominance of the larger magnetic scales, but this remains to 
be demonstrated, and may require proper numerical implementation 
of Braginski viscosity. 

Two types of observations tend to substantiate the necessity for 
this pregalactic mechanism. Beryllium and boron have been detected 
in the atmospheres of very old stars These elements are believed to 
be produced only by spallation of low energy cosmic rays, and thus 
imply the early presence of these cosmic ray^. Further, to produce 
the amount of observed light elements requires cosmic ray intensities 
at this very early time that are comparable to the present intensities. 
This in turn implies the presence of magnetic fields at this early time, 
although, unless most of the interstellar medium at this early time 
was confined to dense clouds, these fields could have been well below 
equipartition with the cosmic rays. 

A second observational result is the detection of Faraday rotation 
by magnetic fields inferred to have microgauss strength, in high-red 
shift-galaxies in an early forming stage. Such galaxies are too young 
for the standard theory to bring seed fields up to such strengths. 

The case for pregalactic mechanisms for magnetic field generation 
is mainly a negative one, based on the nonviability of the standard 
non primordial mechanisms to accomplish the second origin stage. 
The physics of forming galaxies and their discs is not well understood, 
and one has to rely on numerical simulations to draw conclusions 
concerning it. However, these simulations do lead to plausible theories 

^In §2 we mentioned that boron can also be synthesized in supernovae, leaving 
implications of its abundance ambiguous 
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for production of pregalactic magnetic fields by tlie turbulence seen in 
the simulations, provided that the problem of preferred amplification 
at the small scales can be overcome. 

Lastly, the magnetic flux observed in giant radio jets could be dis- 
persed into the intergalactic plasma magnetizing it coherently before 
galaxies form. This provides a very direct pregalactic origin. Unfortu- 
nately, no detailed theory as to how or whether this can happen, has 
yet developed. Also no observational evidence of it has yet emerged. 

In summary, the situation is as follows. The standard theory of 
the generation of magnetic fields through galactic disc dynamos suffers 
from the difficulty of flux expulsion in the case of our Galaxy and no 
consistent theory that resolves this problem has yet been developed. 
A further difficulty is whether the galactic dynamo works on very 
weak fields since most such dynamo theories invoke magnetic fields of 
considerable strength to function. 

Assuming these difficulties cannot be overcome it is necessary that 
there be an initially strong pregalactic field present before the disc 
forms. Such a field can provide the galactic dynamo with a strong 
enough field to function if the fiux expulsion can be handled. 

The flux expulsion might be solved if there are enough superbub- 
bles containing much more than the average number of supernovae 
that actually blow mass entirely out of the galaxy carrying the flux 
with it. The statistics of superbubbles are still unknown and no the- 
ory has yet been developed to show that they can accomplish this 
task. If they can they will produce an inhomogeneous disc fleld that 
is very discontinuous whose lines of force consist of many finite seg- 
ments with ends connected to the intergalactic medium. These are 
the fields that arise if one relies on a galactic wind to be the solution. 
Whether the field arising from such an inhomogeneous theory is com- 
patible with the present field and the astrophysics of the interstellar 
medium associated with it is certainly an open question. 

Because of these serious problems with the disc dynamo it is nec- 
essary to ask whether the pregalactic sources for the disc field are not 
the better solution to the origin problem at least for our Galaxy. Be- 
cause the effort devoted to pregalactic theories is so far inadequate 
to prove that this is a viable origin for pregalactic fields, one cannot 
easily discard this possibility. 

To conclude, we propose the following tentative scenario for the 
origin and growth of galactic magnetic fields is the most likely one: 

1. Fields in the 10~^® - 10~^°G range originate by a battery process 



77 



in cosmological shocks or ionization fronts. 

2. The battery generated fields are sHghtly amphfied by compres- 
sion and significantly amplified by turbulence as protogalaxies 
and similar substructure forms. The resulting fields are in the mi- 
crogauss range, and are coherent over scales of order 1 kpc. This 
requires that amplification be due to the action of the largest 
scale eddies, and that the small scale folded fields produced by 
the small scale eddies quickly saturate. 

3. Once a magnetized galactic disk forms, its magnetic field is sus- 
tained against gas addition and subtraction processes by a su- 
pernova driven a dynamo. Supernova energy and mag- 
netic and cosmic ray buoyancy efficiently expel magnetic flux 
of the wrong sign, allowing the dynamo to work. The field is 
strong enough that amplification at the small scales is no longer 
a problem. 

It appears that to really resolve these questions concerning the 
value of a pregalactic field preceding the for mation of a disc, more 
observations are desperately needed. Hopefully, these can be expected 
in the near future as telescopes and detectors improve. (See Gaensler 
et al. 2004, Gaensler, 2006). They will be able to find many more ex- 
amples of early forming extra galactic systems (damped Lyman alpha 
systems) and determine whether the majority already contain strong 
magnetic fields. 
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